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(^INTRODUCTION 


The  subject  of  our  project  is  diagnostic  imaging  of  breast  tumors.  We  are  testing  hypothesis  that  prevailing 
absorption  of  short  laser  pulses  in  tumors  due  to  advanced  angiogenesis  compared  with  optical  absorption  in 
normal  tissues  will  generate  pressure  profiles  in  breast  tumors  and  convert  these  tumors  into  acoustic  sources. 
Ultrasonic  waves  deliver  diagnostic  information  to  the  breast  surface,  where  signals  are  detected  by  piezoelectric 
transducers.  Reconstruction  of  two-dimensional  images  of  breast  sectors  is  made  on  the  basis  of  time-resolved 
detection  and  analysis  of  wide-band  ultrasound  profiles  recorded  with  a  linear  array  transducer  located  at  the  breast 
surface.  Opto-acoustic  tomography  takes  advantage  of  a  well-developed  signal  processing  and  image 
reconstruction  algorithms,  currently  employed  in  wide-band  ultrasonography. 

The  purpose  of  this  project  is  to  incorporate  laser  excitation  into  a  wide-band  ultrasonography  system  (with 
a  linear  array  probe)  to  yield  a  clinical  prototype  Laser  Optoacoustic  Imaging  System  (LOIS)  and  test  this  system 
initially  in  blood  perfused  phantoms  and  then  in  breast  cancer  patients.  The  tissue  phantoms  will  be  designed  to 
resemble  natural  optical  and  acoustic  properties  of  breast  with  advanced  angiogenesis.  Racial  differences,  such  as 
light  absorption  properties  of  various  skin  types  will  also  be  accounted  for.  Experiments  in  blood  perfused 
phantoms  will  test  the  hypothesis  that  the  breast  circulation  system  and  angiogenesis  can  be  visualized  with  high 
sensitivity  and  contrast.  Clinical  studies  should  confirm  our  hypothesis  that  the  breast  structural  heterogeneity  will 
not  compromise  the  optoacoustic  image  quality,  or  may  be  accounted  for  in  the  optoacoustic  image  processing. 

The  major  anticipated  advantages  of  the  laser  optoacoustic  tomography  are  listed  in  a  sequence 
corresponding  to  their  importance:  (1)  high  optical  contrast  between  tumors  and  normal  tissues  which  results  from 
enhanced  density  of  the  microcirculation  network  within  and  around  tumors  and  enhanced  concentration  of  other 
endogenous  chromophores  in  cancer  cells;  (2)  high  sensitivity  of  detection  which  results  from  efficient  opto- 
thermal  generation  of  pressure  in  tumors,  and  sensitive  detection  of  ultrasonic  waves  with  novel  piezoelectric 
transducers;  (3)  high  (<1  mm)  spatial  resolution  which  results  from  pulsed  laser  excitation  and  time-resolved 
detection  of  undistorted  ultrasonic  wave  profiles  with  array  detectors;  (4)  substantial  depth  of  monitoring  (up  to  8- 
cm)  which  results  from  deep  penetration  of  the  near-infra-red  photons  in  tissues,  efficient  conversion  of  laser 
energy  into  ultrasonic  waves  propagating  through  the  breast  with  insignificant  attenuation  and  minimal  distortion, 
and  from  high  sensitivity  (5  V/bar)  of  piezoelectric  transducers;  and  (5)  applicability  to  dense  breasts. 


(5)  BODY  OF  REPORT 

During  the  Year-1  project  period  all  specific  aims  and  milestones  were  accomplished.  Specific  details  of  the 
optoacoustic  imaging  system  development  and  tests  are  presented  below. 

Task  1  (Year  1):  Development  and  fabrication  of  a  clinical  laser  optoacoustic  imaging  system  (LOIS)  on  the 
basis  of  commercially  available  components  of  a  digital  ultrasound  imaging  system,  such  as  HDI-1000  from  ATL. 


Figure  1. 


Photograph  of  a  clinical  laser  opto-acoustic  imaging  system  LOIS-2. 
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Milestones  of  Task  1  were  accomplished. 

la.  By  the  end  of  the  first  quarter  a  compact  nanosecond  laser  was  incorporated  in  the  optoacoustic 
imaging  system.  The  system  is  made  flexible  to  utilize  either  an  Nd:YAG  laser  (Big  Sky  Lasers,  Inc.  Boseman 
MT)  or  Alexandrite  laser  (CEO  Inc.,  Orlando  FL).  An  HDI-1000  ultrasound  imaging  system  from  ATL  Inc. 
(Seattle,  WA)  was  tested  and  found  not  optimal  as  a  prototype  for  LOIS.  Therefore,  an  electronic  registration 
system  with  low  noise  and  high  dynamic  range  for  fast  signal  acquisition  and  storage  was  developed  by 
LaserSonix  Technologies,  Inc.  (Houston,  TX)  and  provided  to  our  laboratory  as  a  collaborative  effort.  The 
electronic  system  was  then  fabricated  and  assembled  at  the  Pi’s  laboratory  at  UTMB.  By  the  end  of  the  second 
quarter  the  system  was  made  operational  and  tested.  Ultrawideband  piezoelectric  transducer  was  designed  by 
LaserSonix  Technologies,  Inc.  in  collaboration  with  Pi’s  laboratory  and  fabricated  by  Materials  Systems  Inc. 
(Littleton,  MA).  The  electronics  and  software  was  modified  in  accordance  with  the  detection  system  change. 
Software  codes  for  two  dimensional  image  reconstruction  was  completed  and  tested  [1]. 

lb.  By  the  end  of  the  second  quarter  a  software  package  for  automatic  system  control,  signal  processing 
and  image  reconstruction  was  developed  and  tested.  Each  system  component  was  verified  to  possess  optimal 
characteristics  determined  on  the  basis  of  preliminary  results  of  our  experimental  and  theoretical  studies  and  work 
experience  with  existing  bench-top  system.  The  reliability  and  safety  of  the  system  was  tested.  Theoretical 
modeling  of  LOIS  yielded  satisfactory  sensitivity  and  resolution,  so  that  the  clinical  Laser  Optoacoustic  Imaging 
System  (LOIS-2)  was  completed  and  prepared  for  experiments  [2], 

First  experiments  with  LOIS-2  were  performed  in  gel  phantoms  simulating  dense  breast  with  small  tumors 
simulated  by  hemoglobin-colored  spheres.  Two  dimensional  optoacoustic  images  displayed  position  and  and 
diameters  of  spheres  in  agreement  with  their  real  positions  and  diameters  [3,  see  Appendix], 

Further  experiments  with  opto-acoustic  imaging  system  LOIS-2  aimed  at  detection  of  breast  carcinoma 
located  deep  within  a  large  mass  of  breast  tissue.  The  experiment  performed  at  UTMB  demonstrated  potential  of 
LOIS  in  high  contrast  imaging  of  malignant  tumors  [4,5,  see  Appendix].  Breast  specimens  were  obtained  after 
radical  mastectomy  and  examined  within  the  first  ten  minutes  after  surgical  excision.  It  provided  specimens  that 
underwent  minimal  alterations  of  tissue  optical  properties,  except  some  blood  was  drained  form  the  specimens. 
After  opto-acoustic  imaging  procedure,  pathologist  examined  the  mastectomy  specimens  in  order  to  determine 
tumor  geometry.  Location  of  tumors  was  not  known  prior  to  opto-acoustic  imaging  procedure.  Opto-acoustic 
image  of  breast  with  carcinoma  is  presented  in  Fig.  2. 
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Figure  2.  Optoacoustic  image  of  breast  cancer  patient  4081 12P.  Note  that  the  tumor  core  and  sprouts  are 
visualized  with  high  contrast. 

Laser  safety  of  the  system  was  evaluated.  The  opto-acoustic  imaging  procedure  was  performed  using  laser 
fluence  of  50  mJ/cm2  which  corresponds  to  half  maximum  of  permissible  laser  exposure  for  the  skin  in  the  near 
infrared  spectral  range,  as  described  in  American  National  Standards  for  Safe  Use  of  Lasers.  The  laser  irradiance 
of  50  mJ/cm2  can  induce  a  temperature  rise  of  only  50-100  mK  at  the  irradiated  breast  surface  and  about  0.5-1  mK 
in  tumor  located  at  a  5-cm  depth.  The  laser  irradiation  was  performed  with  14-ns  laser  pulses.  Such  a  short 
duration  of  laser  pulses  ensured  that  the  laser-induced  pressure  was  confined  in  the  target  volume  of  tumor  during 
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laser  heating,  yileding  high  resolution  and  and  high  efficiency  of  pressure  generation  with  laser  pulses. 

Optical  fiber  was  scanned  along  the  surface  of  the  radical  mastectomy  specimen  with  4-mm  interval. 
Sixteen  opto-acoustic  signals  were  employed  for  data  collection  and  subsequent  image  reconstruction.  The  scan 
passed  directly  over  the  tumor.  The  tumor  core,  horizontal  tail  and  small  sprouts  may  be  clearly  seen  in  the  opto- 
acoustic  image  (confirmed  by  subsequent  pathology  examination).  The  dimension  of  a  tumor  core  can  be 
estimated  as  10x7  mm.  Tumor  was  detected  at  the  depth  of  1 1  mm  from  the  surface.  Skin  was  removed  from  most 
part  of  the  specimen,  therefore  the  irradiated  surface  is  depicted  as  a  rough  line.  Opto-acoustic  contrast  between 
the  tumor  and  surrounding  normal  tissues  in  this  image  exceeds  350%  making  surrounding  normal  tissue  looking 
dark  on  the  image. 

The  same  breast  was  imaged  in  vivo  with  General  Electric  mammography  system  with  Mo  anode  and  Mo 
filter  and  with  Acuson  128  XP/10  system  employing  3.5  MHz  C3-curved  linear  array  transducer.  These  images 
are  depicted  in  Fig.  3.  X-ray  image  revealed  all  details  of  this  tumor  including  core  with  microcalcifications  and 
lateral  sprouts.  However,  x-ray  also  showed  other  areas  of  the  breast  similarly  bright  on  mammography  image, 
however,  found  normal.  The  ultrasound  image  of  the  same  tumor  was  not  very  impressive.  The  core  of  tumor  and 
its  long  tail  may  also  be  seen  in  this  image.  However,  small  sprouts  are  not  visible  to  limited  acoustic  contrast  [5]. 


Figure  3.  X-ray  radiography  image  obtained  in  vivo  with  General  Eelectric  mamography  system.  Mo  anode, 
Mo  filter  (right)  and  ultrasound  image  obtained  with  Acuson  128  XP/10  3.5  MHz  USI  system 
(left)  of  a  malignant  tumor  in  the  breast.  UTMB  Patient  #  4081 12P.  Note  that  contrast  of  the  tumor 
relative  to  background  is  not  sufficient  to  define  exact  shape  and  dimensions  of  this  malignancy. 

Using  the  experimental  data  obtained,  one  can  calculate  minimal  radius  of  spherical  tumors  with  various 
blood  contents,  which  can  currently  be  detected  inside  the  breast  [5-7].  The  radius  of  the  minimal  detectable 
spherical  tumor  is  plotted  in  Fig.  4  as  a  function  of  depth  from  the  irradiated  surface  for  various  absorption 
coefficients  corresponding  to  different  blood  content  in  tumors.  The  anticipated  blood  content  in  actively  growing 
small  tumors  is  about  4-16%  as  found  in  tumor  pathology  studies. 


Figure  4.  Minimally  detectable  tumor  radius  as  a  function  of  tumor  depth  in  the  bulk  of  normal  tissue. 
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v  The  areas  above  the  curves  in  Fig.  4  represent  parameters  of  tumors,  which  can  be  detected  by  wide-band 
acoustic  transducers  made  of  piezoelectric  composite  materials  used  in  our  experiments.  The  areas  under  the 
curves  represent  areas  below  the  limit  of  detection  for  these  transducers.  Fig.  4  demonstrates  that  the  sensitivity  of 
LOIS-2  permits  detection  of  2-mm  tumors  with  an  average  blood  content  of  8%  at  a  depth  of  ~6  cm.  Therefore, 
the  sensitivity  of  our  current  laser  optoacoustic  imaging  system  is  sufficient  for  detecting  small  tumors  (with  or 
without  calcification)  using  a  safe  level  of  laser  irradiation.  The  limits  of  detection  for  laser  optoacoustic 
tomography  may  be  further  improved  with  the  use  of  state-of-the-art  acoustic  transducer  arrays  supported  by  the 
most  advanced  electronics  and  image  processing  software.  This  indicates  that  LOIS  has  a  potential  for  early 
detection  of  breast  cancer,  which  will  be  tested  in  later  stages  of  our  project. 


Task  2  (Year  1):  Evaluation  of  the  laser  optoacoustic  imaging  system  using  tissue  phantoms  simulating  dense 
breast  with  advanced  tumor  angiogenesis. 

Milestones  of  Task  2  were  accomplished. 

Milestone  2a.  By  the  end  of  the  third  quarter  the  clinical  system  performance  was  tested  in  tissue 
phantoms  resembling  optical  and  mechanical  properties  of  breast  with  tumor  angiogenesis.  The  clinical  system 
sensitivity,  the  accuracy  of  the  “tumor”  localization,  the  image  contrast  and  spatial  resolution  will  be  compared 
with  the  state-of-the-art  ultrasound  and  mammography  systems. 

Clinical  optoacoustic  imaging  system  was  assembled  and  its  performance  was  evaluated  in  comparison  with 
our  current  bench-top  system  using  tissue  phantoms  made  of  chicken  breast  slabs  and  milk-gelatin  slabs  with 
embedded  blood  vessels  extracted  from  mice  and  rats  and  operated  using  closed  loop  peristaltic  pump  [4,8,9], 

Figurre  5  (a,b)  shows  two  optoacoustic  images  of  a  gel  phantom  with  optical  properties  resembling  dense 
breast  with  one  artery  and  one  vein.  Blood  vessels  surgically  extracted  from  a  rat  and  filled  with  blood,  one  with 
deoxygenated  blood  and  one  with  oxygenated  blood,  we  embedded  inside  the  gel  phantom.  A  cross-sectional  two- 
dimensional  image  was  made  of  this  tissue  phantom.  External  diameter  of  blood  vessels  was  approximately  0.72 
mm  and  the  lumen  diameter  was  0.44  mm. 
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Figure  5.  Optoacoustic  images  of  tissue  phantoms  made  of  milk  gelatin  with  embedded  two  blood  vessels 
filled  with  blood.  Imaging  was  performed  with  two  different  laser  wavelengths:  532  nm  (a)  and 
1064  nm  (b).  One  blood  vessel  had  oxygen  saturation  level  of  hemoglobin  equal  70%.  The  other 
blood  vessel  had  oxygen  saturation  level  of  about  10%. 

One  can  see  that  images  of  both  blood  vessels  have  equal  brightness  at  the  wavelength  of  532  nm  where 
oxyhemoglobin  and  deoxyhemoglobin  have  similar  optical  absorption  (Fig.  5a),  the  effect  anticipated  also  at  the 
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wavelength  of  Alexandrite  laser  operating  at  780  nm.  However,  at  the  wavelength  of  1064  nm  blood  vessel  with 
oxygenated  arterial  blood  has  much  brighter  image  compared  with  the  blood  vessel  filled  with  deoxygenated  blood 
(almost  invivible)  in  Fig.  5b  [4,7,8]. 

The  same  phantom  was  imaged  with  General  Electric  mammography  system  anode  and  filters  made  of 
molibdenum  and  with  Acuson  128  XP/10  system  employing  3.5  MHz  C3-curved  linear  array  transducer.  This 
experiment  was  designed  and  performed  in  order  to  compare  the  optoacoustic  system  sensitivity,  the  accuracy  of 
the  blood  vessel  localization,  the  image  contrast  and  spatial  resolution  with  those  of  the  state-of-the-art 
mammography  and  ultrasonography  systems.  However,  neither  x-ray  mammography  nor  ultrasonography  could 
visualize  blood  vessels  inside  dense  gelatin  phantom.  The  reason  for  X-ray  system  failure  is  well  known  and  can 
be  explained  by  high  radiologic  density  of  phantom  tissue  (similar  to  density  of  a  young  woman).  The  reason  for 
ultrasound  failure  is  the  absence  of  acoustic  contrast  in  density  of  blood  vessels  and  surrounding  gelatin.  This 
experiment  reveal  potential  advantage  of  optoacoustic  imaging  system  compared  with  existing  modalities. 

Milestone  2b.  By  the  end  of  the  fourth  quarter  the  optimal  characteristics  of  each  system  component  from 
the  stand  point  of  diagnostic  efficacy,  and  safety  was  determined.  The  system  parameters  were  improved  utilizing 
the  results  of  the  preliminary  evaluation  in  tissue  phantoms.  The  changes  in  the  transducer  array  design  suggested 
by  results  of  in  vitro  studies  were  incorporated. 


Figure  4. 


Photograph  of  the  linear  array  of  32  ulrawideband  ultrasonic  transducers  with  improved 
characteristics:  sensitivity  5  V/bar,  resolution  0.5x1  mm  in  wide  focal  zone. 


In  order  to  optimize  optoacoustic  system  characteristics  form  the  stand  point  of  diagnostic  efficacy,  safety, 
resolution  and  sensitivity,  such  laser  irradiation  parameters  as  dynamic  range  and  ultrasonic  bandwidth  of  the 
acoustic  detection  were  varied.  Optimal  parameters  of  the  array  transducer  were  found  to  be  as  follows:  (1) 
ultrasound  detection  bandwidth:  20  kHz  to  4  MHz,  dimensions  of  each  transducer  1.5  mm  x  1mm  with  1  mm  gaps, 
6  transducers  in  one  piezoelement,  64  piezoelements  total  in  120°  linear  array  with  diameter  of  120  mm.  Sensitivity 
of  such  transducer  made  of  PZT-5H  ceramic  embedded  in  polymer-epoxy  resin  with  2-2  structure  was 
exceptionally  high,  5  V/bar  with  effective  noise  pressure  of  only  4  pbar  [6,10]. 


(6)  KEY  ACCOMPLISHMENTS 

•  Laser  Opto-acoustic  imaging  system  (LOIS)  was  developed  and  assembled  for  two-dimensional  tomography  of 
early  breast  cancer  utilizing  components  of  wide-band  ultrasound  imaging  system  and  nanosecond  Nd:YAG  or 
Alexandrite  lasers. 

•  LOIS  was  tested  in  gelatin  phantoms  simulating  dense  breast  with  blood  circulation.  Experiments  demonstrated 
that  LOIS  is  capable  of  not  only  visualizing  blood  vessels,  but  also  differentiating  veins  with  deoxygenated  blood 
from  arteries  with  oxygenated  blood. 

•  LOIS  was  tested  in  real  breast  with  carcinoma  ex- vivo  (radical  mastectomy  specimen).  Experiment  demonstrated 
that  LOIS  can  visualize  details  of  malignant  tumor  not  present  in  ultrasound  images. 
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(8)  CONLUSIONS 

•  Laser  opto-acoustic  imaging  system  combining  pulsed  optical  excitation  of  pressure  profiles  and  sensitive  ultrawide- 
band  detection  of  ultrasound  is  capable  of  visualizing  breast  tumors  with  exceptional  contrast  relative  to  normal  breast 
tissue  background. 

•  An  evidence  based  on  imaging  blood  vessels  with  oxy-  and  deoxy-  blood  was  obtained  that  LOIS  may  be  capable  of 
differentiating  malignant  tumors  (with  deoxy-blood  in  tumor  angiogenesis)  from  benign  tumors  (with  oxy-  blood  in 
tumor  angiogenesis). 

•  Optimal  parameters  for  clinical  LOIS  are  as  follows: 

Optical  Excitation  -  10  ns  to  100  ns  laser  pulses  with  repetition  rate  of  10  Hz,  fluence  of  50  mJ/cm2, 
and  wavelength  chosen  from  the  near  infrared  range  between  680  nm  and  1064  nm  and  dependent  on  the  type  of 
tumor  to  be  imaged. 

Ultrasonic  Detection  -  piezoelectric  transducers  made  of  2-2  ceramic-polymer  composite  material 
operating  in  an  ultrawide  band  of  ultrasonic  frequencies  from  20  kHz  to  4  MHz. 
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ABSTRACT 


The  second  generation  of  the  laser  optoacoustic  imaging  system  (LOIS-2)  for  breast  cancer  detection,  localization  and 
characterization  using  a  32-element  arc-shaped  transducer  array  was  developed  and  tested.  Each  acoustic  transducer  was 
made  of  llO-pm  thick  SOLEF®  PVDF  film  with  dimensions  of  1mm  x  12.5  mm.  The  frequency  band  of  transducer  array 
provided  0.4-mm  axial  in-depth  resolution.  Cylindrical  shape  of  this  10-cm  long  transducer  array  provided  an  improved 
lateral  resolution  of  1.0  mm.  Original  and  compact  design  of  low  noise  preamplifiers  and  wide  band  amplifiers  was 
employed.  The  system  sensitivity  was  optimized  by  choosing  limited  bandwidth  of  ultrasonic  detection  20-kHz  to  2-MHz. 
Signal  processing  was  significantly  improved  and  optimized  resulting  in  reduced  data  collection  time  of  13  sec.  The 
computer  code  for  digital  signal  processing  employed  auto-gain  control,  high-pass  filtering  and  denoising.  An  automatic 
recognition  of  the  opto-acoustic  signal  detected  from  the  irradiated  surface  was  implemented  in  order  to  visualize  the  breast 
surface  and  improve  the  accuracy  of  tumor  localization.  Radial  back-projection  algorithm  was  used  for  image  reconstruction. 
Optimal  filtering  of  image  was  employed  to  reduce  low  and  high  frequency  noise.  The  advantages  and  limitations  of  various 
contrast-enhancing  filters  applied  to  the  entire  image  matrix  were  studied  and  discussed.  Time  necessary  for  image 
reconstruction  was  reduced  to  32  sec.  The  system  performance  was  evaluated  initially  via  acquisition  of  two-dimensional 
opto-acoustic  images  of  small  absorbing  spheres  in  breast-tissue-like  phantoms.  Clinical  ex-vivo  studies  of  mastectomy 
specimen  were  also  performed  and  compared  with  x-ray  radiography  and  ultrasound. 

Keywords:  Opto-acoustic  imaging,  breast  cancer,  transducer,  laser 


1.  INTRODUCTION 

In  previous  publications14  we  described  the  basic  ideas,  principles  and  the  first  clinical  prototype  of  the  laser 
optoacoustic  imaging  system  (LOIS-01).  This  system  LOIS-01  was  assembled  of  a  compact  Nd:YAG  laser,  fiberoptic  light 
delivery  system,  a  linear  array  of  12  wide-band  acoustic  transducers,  and  a  data  acquisition  card  operated  by  computer. 
LOIS-01  was  successfully  tested  in  imaging  mastectomy  specimen.  Total  of  fifteen  radical  mastectomy  specimens  was 
investigated  at  UTMB  and  UT/MDACC.  In  all  of  15  cases  the  tumor  location  was  detected  and  its  depth  and  relative 
position  were  measured.  Comparison  of  opto-acoustic  images  with  those  of  X-ray  radiography,  MRI  and  ultrasound  imaging 
performed  in  vivo  revealed  good  correlation  in  tumor  size  and  location.  A  pronounced  opto-acoustic  contrast  of  200-450% 
between  breast  tumors  and  normal  breast  tissues  was  demonstrated.  This  contrast  substantially  exceeds  any  other 
endogenous  tissue  contrast  currently  utilized  in  clinical  ultrasonography,  MRI  and  X-ray  mammography.  Based  on  literature 
data  and  our  gross  observations  of  tumor  cross-sections  we  hypothesize  that  the  opto-acoustic  contrast  results  primarily  from 
increased  optical  absorption  in  the  dense  micro  vascularity  of  the  tumors  . 
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In  patients  receiving  radiotherapy,  tumors  were  found  to  contain  enhanced  concentration  of  dense  highly  scattering 
fibrotic  tissue.  Especially  designed  for  in  vitro  studies  LOIS-01  was  characterized  by  relatively  high  sensitivity  (30  |TV/Pa), 
but  was  rather  cumbrous.  A  small  number  of  acoustic  transducers  (12)  in  receiving  linear  array  limited  the  lateral  resolution 
of  images.  The  objective  of  the  present  study  was  to  fabricate,  assemble  and  test  the  prototype  of  laser  opto-acoustic  imaging 
system  (LOIS-02)  for  clinical  investigations.  Design  of  the  new  system  had  to  be  compact  and  suitable  for  performing 
studies  in  vivo.  The  resolution  of  acoustic  transducer  array  had  to  be  improved  for  small  tumor  detection  and  measurement  of 
its  dimensions.  New  signal  and  image  processing  employed  had  to  allow  operator  to  reconstruct  tumor  image  in  a  short  time 
after  data  collection. 


2.  MATERIALS  AND  METHODS 


2.1.  Arc  array  transducer 

A  specially  designed  array  of  PVDF  acoustic  transducers  was  employed  in  LOIS-02.  The  array  had  32  rectangular 
piezoelectric  transducers  of  1x12.5  mm  dimensions  and  3.85  mm  distance  between  transducers  (see  Fig.l).  Piezoelectric 
polymer  PVDF  of  llO-pm  thickness  was  used  for  transducer  fabrication.  Low  acoustic  impedance  and  ability  to  operate  in 
wide  ultrasonic  frequency  band  are  the  advantages  of  PVDF  for  detection  of  opto-acoustic  profiles.  The  transducers  were 
located  on  arc  surface  of  60-mm  radius.  This  geometry  provided  a  high  resolution  in  the  vicinity  of  focal  zone.  Backing 
material  was  used  for  mechanical  matching  of  transducer  and  damping  reverberations  after  detection.  When  transducer  is 
excited  by  acoustic  pulse  it  vibrates  at  its  natural  resonance  frequency. 


Figure  1.  Schematic  diagram  of  arc  transducer  array. 


If  acoustic  impedance  of  backing  material  is  similar  to  that  of  transducer  the  ringing  could  be  damped  out  and  the 
frequency  bandwidth  increased.  The  front  surface  of  transducer  array  was  covered  by  aluminum  film  to  reduce  external 
electric  noise.  The  thickness  of  Aluminum  layer  was  optimized  for  the  best  relation  of  ultrasound  penetration  and  electric 
noise  reduction.  Transducer  array  housing  also  was  made  of  aluminum.  Acoustic  isolating  material  was  placed  between 
piezoelectric  elements  and  housing  to  prevent  ringing  of  the  housing  that  follows  the  vibration  of  the  piezoelectric  elements 
themselves. 

Absolute  sensitivity  of  each  transducer  was  measured.  Results  of  this  calibration  are  presented  in  Fig.  2.  The  mean 
sensitivity  value  for  transducer  array  of  5  =  6  pV/Pa  is  three  times  lower  than  the  ultimate  theoretical  value.  The  ultimate 
theoretical  sensitivity  of  transducers  can  be  calculated  on  the  basis  of  piezoelectric  properties  of  PVDF  material  employed: 


$  ~  833d  • 


(1) 


37 


where  g33  is  the  piezoelectric  pressure  constant  of  material,  d  is  the  thickness  of  PVDF  film. 
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Figure  2.  The  sensitivity  of  transducers  in  the  array.  Maximum  corresponds  to  the  sensitivity  of  7  ^V/Pa. 
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of  capacitance  of  piezoelectric  element.  The  capacitance  CT  can  be  found  as: 

cy  =  Hsd 

Ti  (2) 

iTlow  toTZ  "“I  °f  8  Pi“°eTM‘  “•  “  d'C,Cmk  °“r  of  each  transducer 

P  “ar^rf  e  £  I  ,  ' ‘  “  ""  ”*y  Capad,a"Ce  °f  COnn“»ng  Wires  *”d  °p  el=«o»ic 

preamphfier.  The  shunt, „g  effect  ,s  the  source  of  additional  loses  of  sensitivity  in  the  system.  I,  is  difficult  calculate  all 
possible  losses  of  acousttc  energy,  therefore  the  direct  measurements  of  transducers  sensitivity  in  cay  is  of  great  importance 

The  vartabtitty  on  the  sensitivity  between  array  elements  is  about  15*  due  to  the  lack  of  reproducibility  in  the  manual 
assembly  process. 

2.2.  Sensitivity  of  small  tumor  detection 

One  of  the  important  parameters  of  LOIS  is  the  ultimate  sensitivity  for  small  tumor  detection.  Let  us  consider  a 

umor  as  a  sphere  of  radius  o  with  optical  absorption  coefficient  Ma.  The  maximum  depth  of  the  tumor  detection  can  be 
evaluated  as! 


=  -— Inf  -iM- 
Veff  I>aaF0 


where  fieff  is  the  effective  light  attenuation  coefficient  in  breast  tissue,  r  is  the  distance  between  the  tumor  and  transducer. 

laser  proems  ^  T  ^  mUSt  *  ^  ***  01  ***  due  t0  safety  requirement  for  medical 

P  '  e  Parameter  r  characterizes  the  efficiency  of  thermo-acoustic  excitation  in  the  media.  For  live 

biological  tissues  at  37  °C  this  parameter  varies  in  the  range  0.40  -  0.85  being  maximal  in  fatty  tissues.  The  effective  noise 
pressure  pN  »  calculated  for  signal-to-noise  ratio  (SNR)  equal  1.  The  root-mean-square  of  noise  voltage  URMS  was  measured 
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for  each  transducer.  This  value  is  defined  by  the  thermal  noise  of  transducer  capacitance  and  input  noise  of  preamplifiers. 
For  LOIS-2  the  measured  value  of  URMS  was  about  30  jiV.  The  noise  voltage  is  reduced  "JIN  times  if  the  procedure  of 
averaging  by  N  pulses  is  employed.  We  used  N=16  in  our  measurements  that  allowed  us  to  get  value  of  effective  noise 
pressure  as  1.3  Pa.  Substituting  pN  in  formula  (3)  yields  the  maximum  depth  of  tumor  detection  of  about  6  cm  assuming  the 
tumor  to  be  a  sphere  of  2  mm  diameter  with  pa  =1  cm*1. 


2.3.  Lateral  resolution 

A  spatial  distribution  of  acoustic  field  in  the  focal  surface  of  array  can  be  described  by  corresponding  expression  for 
the  focused  arc  transducer3: 


Pf  =  P/0 


kxQ 


max 


(4) 


where  pp  is  the  pressure  amplitude  in  the  focal  zone,  k  is  the  wave  number,  ©j^  =  arcsin(r/P)  is  the  angular  aperture 

of  array,  defined  by  the  radius  of  arc  curvature  R  and  the  aperture  T  (see  Fig.4).  The  size  of  focal  waist  <5/  defines  the  lateral 
resolution  of  the  array.  A  value  of  focal  waist  can  be  found  from  the  relation  (4)  as  the  distance  between  zeroes  of  pressure 
of  -1  and  +1  orders: 


1.22 


A 

arcsin(r/i?) 
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where  X  is  the  wavelength  of  ultrasound  wave. 


Figure  4.  Schematic  diagram  of  directivity  of  acoustic  detection  for  an  arc-shaped  transducer. 

For  our  transducer  array  the  ratio  r/Re quals  approximately  0.87,  therefore  the  expression  (5)  can  be  simplified, 
yielding  <5/  =  1.22A.  This  estimation  is  valid  for  the  monofrequency  acoustic  waves.  In  case  of  pulsed  detection  it  is 
necessary  to  take  into  account  the  frequency  spectrum  and  to  choose  the  value  of  effective  X  correctly.  The  expression  (5) 
provides  correct  results  when  all  transducers  in  array  are  involved  in  focused  pattern  formation.  A  rectangular  transducer 
with  side  dimension  b  has  directivity  pattern  with  angular  width  (X  measured  at  -3  dB  level  given  by 

(X  =  arcsin(0.45A/6)  (6=1  mm  for  the  transducer  employed  in  our  studies).  This  means  that  the  lower  the  ratio  of  X/b 
(or  the  lower  the  value  of  A  with  fixed  value  of  6),  the  narrower  the  directivity  pattern  of  a  single  transducer  and  smaller  the 
space  with  resolution  defined  by  the  value  of  A.  Therefore  the  resolution  depends  on  the  position  of  tumor  relative  to  the 
array.  In  the  vicinity  of  the  focal  zone  (i.e.  for  the  points  of  location  inside  the  circle  of  1-cm  radius  with  center  in  the 
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geometrical  focus  of  the  array)  the  resolution  is  less  1-mm.  In  the  main  part  of  space  the  resolution  is  about  Imm<<5/<2  mm. 


2.4.  Axial  in  depth  resolution 

The  axial  (in-depth)  resolution  of  the  pulsed  opto-acoustic  system  operating  under  irradiation  conditions  of  temporal 
pressure  confinement  is  directly  defined  by  the  detection  bandwidth  of  acoustic  transducers.  Two  coherently  emitting 
spherical  acoustic  sources  with  radius  a  will  be  detected  as  separate  objects  if  the  space  between  their  centers  equals  or 
greater  than  Az  =  la.  For  correct  detection  of  pulsed  opto-acoustic  waveform  excited  in  absorbing  sphere  of  radius  a,  it  is 
necessary  to  detect  opto-acoustic  signals  within  a  wide  ultrasound  frequency  band  of  0  -fb1: 


fb  =  0.15cs/  a  (6) 

where  Cs  is  the  speed  of  sound,  fb,  is  the  upper  limit  of  the  ultrasonic  frequency  band  in  the  frequency  spectrum  of  N-shaped 
signal  emitted  by  the  small  sphere. 

Therefore  the  following  inequality  must  be  satisfied:  fb  <  where  fmax  is  the  upper  frequency  limit  of 


ultrasound  detector.  Using  expression  (6)  for  fb ,  one  can  obtain: 


Az  >  1.5 


fm 


(7) 


The  upper  frequency  limit  of  ultrasonic  detection  band  fmax  is  defined  by  the  transducer  thickness  and  the  frequency 
bandwidth  of  preamplifier  employed.  In  our  design  the  upper  frequency  limit  was  6  MHz,  that  provided  0.38  mm  axial  in- 
depth  resolution. 


2.5.  LOIS-02  Description 

Schematic  diagram  of  LOIS-02  is  shown  in  Fig.  5.  A  Nd:YAG  laser  (Big  Sky  Lasers,  MA)  operating  at  the 
wavelength  of  1064  nm  was  used  as  a  source  of  near-infrared  pulses  of  10-ns  duration.  A  quartz  optical  fiber  was  employed 
for  the  laser  pulse  delivery  to  the  tissue  surface.  The  fiber-optic  system  was  attached  to  a  three-dimensional  translation  stage 
driven  by  a  step  motor.  An  electronic  data  acquisition  system  was  synchronized  with  the  linear  translation  along  the  chosen 
direction. 


Figure  5.  Schematic  diagram  of  LOIS-02. 


Breast  phantoms  and  mastectomy  specimens  were  placed  on  the  surface  of  the  acoustic  array  with  m=32  elements. 
Transducer  array  housing  was  fixed  vertically  in  a  specially  designed  holder,  which  provided  a  convenient  access  to  the 
samples.  A  circuit  board  with  32  identical  preamplifiers  was  placed  in  the  same  housing  in  the  vicinity  of  transducers.  Each 
preamplifier  had  a  1  MOhm  electrical  input  impedance  and  low  output  impedance  of  50  Ohm.  It  allows  one  to  match  the 
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output  impedance  of  piezoelectric  elements  and  the  connecting  cable  in  order  to  improve  the  sensitivity  and  the  signal-to- 
noise  ratio.  The  voltage  of  electrical  signals  from  preamplifiers  was  amplified,  and  signals  were  then  digitized  with  a  PDA- 
12  data  acquisition  card  with  two  input  channels  (SignaTech,  CA).  A  multiplexor  sequentially  connected  thirty-two  amplifier 
outputs  with  the  two  input  channels  of  the  data-card.  Acquired  data  were  stored  in  a  computer  for  further  processing. 

2.6.  Data  processing  in  LOIS-02. 

Signal  processing  algorithm  realized  in  LOIS-02  is  shown  schematically  in  Fig.6.  Each  signal  detected  by  acoustic 
transducers  in  the  array  consists  of  the  following  three  components:  (1)  a  sharp  peak  produced  at  the  very  surface  of  the 
medium  under  investigation,  (2)  a  low-frequency  smooth-rise  exponential  slope  produced  by  the  attenuation  of  light  intensity 
inside  tissue,  and  (3)  short  N-shaped  pulses  generated  by  small  absorbing  tumors.  We  assumed  here  that  the  tumor  is  a 
sphere  of  radius  a .  The  duration  of  the  N-shaped  pulse  is  defined  by  the  time  of  sound  propagation  through  the  sphere 
diameter.  The  delay  in  time  of  signal  arrival  depends  on  tumor  location  relative  to  transducer.  Therefore  position  of  the 
tumor  and  its  dimensions  can  be  determined  from  the  opto-acoustic  signal. 


Before  the  signal  could  be  used  for  image  reconstruction,  it  is  necessary  to  eliminate  completely  its  low  frequency 
component  associated  with  homogeneous  light  attenuation  inside  the  tissue,  because  it  can  significantly  decrease  the  image 
contrast.  This  was  achieved  in  LOIS  with  high-pass  numeric  filtering  with  RC  cells.  The  value  of  RC-10  ps  defines  the  cut¬ 
off  ultrasonic  frequency,  and  the  sharpness  of  the  filter  cut-off  is  defined  by  the  number  of  RC  cells.  The  greater  the  number 
of  elementary  RC  cells,  the  steeper  the  slope  of  the  filter  transfer  function.  The  cut-off  frequency  and  the  slope  of  the  transfer 
function  of  the  filter  could  be  varied  conveniently  by  the  imaging  operator. 

The  signal  generated  by  laser  pulses  in  the  subsurface  layer  (very  close  to  the  surface)  possesses  maximum 
amplitude  over  the  entire  detected  course.  Therefore,  the  signal  gradient  is  maximal  at  the  tissue  surface,  which  yields  high 
ultrasound  frequencies.  High  ultrasound  frequencies  that  fall  near  the  resonance  frequency,  /max,  of  piezo-elements  can 
induce  significant  reverberation  in  the  transducer.  Our  transducers  were  designed  to  damp  resonance  frequencies  and  widen 
the  range  of  detectable  frequencies.  Nevertheless,  a  strong  reverberation  takes  place  after  the  arrival  of  the  signal  from  the 
surface  of  tissue.  These  reverberations  do  not  affect  the  useful  part  of  the  opto-acoustic  profile.  Besides,  it  is  convenient  to 
employ  the  sharp  peak  with  high  gradient  (derivative)  for  automatic  determination  of  the  tissue  surface.  The  position  of 
absolute  minimum  in  the  signal  was  considered  as  the  position  of  the  irradiated  tissue  surface.  The  controlling  computer 
generated  numerical  signal  at  the  surface  position  determined  for  each  transducer  and  stored  these  data  in  a  separate  file.  All 
data  samples  were  cut  off  the  signal  after  arrival  of  the  signal  generated  at  the  tissue  surface.  Integrated  signals  were  used  as 
input  data  for  the  image  reconstruction  code.  The  surface-marks  generated  by  computer  were  used  for  the  surface 
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visualization.  Combined  image  contained  the  surface  and  the  tumors,  provided  information  about  tumor  location  relative  to 
the  irradiated  surface. 


2.7.  Image  reconstruction  algorithm 

Let  us  consider  the  case  of  thermo-optical  excitation  of  acoustic  waves  in  homogeneous  medium  with  coefficient  of 
optical  absorption  fiQ,  The  temporal  integral  of  acoustic  pressure  detected  by  transducer  located  in  the  point  r  can  be 
expressed  in  the  following  form6: 

I  r  —  r'  | 

,  a  Hatrmnut-i — b 

u(r,t)  =  j  p(r,t')dt'  =  — — J - i"z  ~ , - — - d? .  (8) 

L  4nCpi  \r~r\ 

where  (5  is  the  thermal  coefficient  of  volume  expansion,  Cp  is  the  specific  thermal  capacity  of  the  medium,  L  is  the  laser 
waveform  function.  All  other  notations  used  in  formula  (8)  are  shown  in  Fig.7. 


Figure  7.  Geometry  of  the  opto-acoustic  detection. 


Integral  is  calculated  over  the  entire  space.  It  means  that  the  acoustic  pressure  p(r,t)  at  the  time  moment  t  and  in 
the  point  r  is  determined  by  the  thermal  sources  located  in  the  spherical  layered  volume  with  radius  |  T  —  r'  |  and  thickness 
dr' .  Acoustic  waves  arrive  to  the  measurement  point  r  with  time  delay  |  r  —  7  \  jcs  .  If  the  laser  pulses  are  sufficiently 
short,  then  generation  of  thermal  sources  may  be  considered  instantaneous  and  the  laser  waveform  function  can  be  expressed 
in  the  form:  L(t)  =  TLS(t) ,  where  rL  is  the  laser  pulse  duration,  8(t)  is  the  delta  function.  Equation  (8)  for  short  laser  pulses 
can  be  written  in  the  following  form: 


u(r,t) 


fi  (r'  )F(r'  )S(t  -  — — — ) 

P  f _ 

4  nCp\  jr-F'| 


dr' 


(9) 


where  F  is  the  laser  fluence.  Let  the  surface  2  be  the  spherical  surface  of  radius  R  with  origin  at  the  position  of  receiving 
transducer.  Radius  of  the  sphere  is  defined  by  the  time  of  sound  propagation  from  the  surface  to  the  transducer.  Vector  R 
can  be  presented  as  R  =  Cstn  ,  where  n  is  the  unit  vector  of  the  surface  element  dL .  Thus,  the  temporal  integral  of  acoustic 
pressure  can  be  expressed  as: 

u(r,t)  =  —@- - f  Q(r  +  cstn)d?:,  (10) 

4  nCpcjt  Jz 

where  function  Q(r)  =  fla(r)F(r)  describes  the  spatial  distribution  of  thermo-acoustic  sources. 
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Equation  (10)  can  be  considered  as  the  fundamental  equation  for  the  opto-acoustic  tomography.  The  temporal 
integral  of  pressure  profile  detected  by  transducer  at  the  time  t  is  the  superposition  of  thermal  sources  located  on  the  sphere  of 
radius  R  =  Cst.  The  product  of  light  absorption  coefficient  and  light  fluence  at  this  spherical  surface  determines  acoustic 

amplitude  of  the  thermal  source. 

2.8.  Image  processing 

An  image  reconstruction  algorithm  based  on  Radon  transform  with  limited  set  of  data  was  proposed  for  the  opto- 
acoustic  tomography7.  However,  this  method  does  not  yield  satisfactory  resolution  even  with  significant  number  of 
transducers  employed.  We  adopted  a  radial  back-projection  algorithm,  earlier  developed  for  the  ultrasound  imaging8.  For  the 
reconstruction  of  two-dimensional  opto-acoustic  images  we  used  signal  integrals  shown  in  Fig  6  and  projected  them  back 
onto  the  two-dimensional  grid  taking  into  account  the  directivity  pattern  (angle  of  acceptance)  of  each  transducer.  The  entire 
image  field  resulted  from  radial  back-projection  algorithm  can  be  written  in  the  following  form. 


171  ,  .  .  , 

£(r)=  lun{\r-rn\/cs)\r-7n\  OD 

n= 1 

where  un  is  the  integrated  signal  of  the  n-th  transducer,  rn  is  the  radius  vector  of  the  n-th  transducer,  m  is  the  number  of 
transducers  in  array.  The  image  represents  distribution  of  a  product  of  thermo-acoustic  efficiency,  optical  absorption  and 
absorbed  laser  energy.  The  back-projected  images  are  noisy  and  their  contrast  is  insufficient  for  correct  object  recognition. 
Therefore  an  image  filtration  procedure  was  employed: 

^(r)  =  J—\FA(b)H{(b)emda)  (12) 

4  n  J  s 

where  F^{(b)  =  J  £(r)e  d(b  is  the  2D  Fourir  spectrum  of  the  back-projected  image,  H((b)  =  \a>\e  ^  is  the 

filter  transfer  function,  (b  =  {  0)x , Q)y  }  is  the  vector  of  spatial  frequencies.  The  parameter  <7  allows  us  to  modify  the  filter 

bandwidth  and  the  steepness  of  filter  characteristic.  It  can  be  done  manually  by  operator  or  automatically.  In  automatic 
procedure  the  parameter  C  was  chosen  if  additional  condition  for  image  quality  was  applied.  The  filter  parameter  O  was 
optimized  employing  the  criterion  of  maximum  contrast  of  the  resulting  image. 


2.9.  Time  parameters  of  LOIS-02. 

Data  acquisition  time  in  LOIS-02  equals  13  sec  for  one  fiber  position.  The  laser  pulse  repetition  rate  is  assumed  to 
be  20  Hz  and  16  pulses  for  averaged  acquisitions  are  employed.  Four  seconds  is  needed  for  data  processing  and  back- 
projected  image  formation.  Image  filtration  and  its  optimal  processing  takes  10  sec.  Total  time  equals  32  sec. 


3.  RESULTS  AND  DISCUSSION 


3.1.  Spheres  in  milk 

Laser  optoacoustic  imaging  system  was  tested  in  various  phantoms.  Two-dimensional  image  of  a  sphere  8-mm  in 
diameter  in  milk  is  shown  in  Figure  8.  Whole  milk  was  diluted  by  water  to  obtain  the  effective  optical  attenuation  coefficient 
1.2  cm'1  to  mimic  the  optical  properties  of  the  breast  tissue.  Optical  absorption  coefficient  of  the  sphere  was  about  1.0  cm’ . 
This  value  is  typical  for  the  light  absorption  in  tumors  we  measured  in  mastectomy  specimen4.  The  sphere  was  located  near 
the  focal  zone  of  transducer  array.  Axes  in  images  are  scaled  in  millimeters.  Position  of  the  transducer  array  is  marked  by 
arc  dotted  line.  Surface  position  was  detected  and  reconstructed  by  using  the  computer-generated  surface  marks.  In  this 
particular  case  the  surface  was  flat,  therefore  it  was  marked  by  a  straight  line.  A  real  surface  position  (72-mm)  was  m 
agreement  with  automatically  detected  value.  The  center  sphere  depth  (19-mm)  in  image  corresponded  very  well  with 
measured  value. 


43 


Figure  8.  Opto-acoustic  image  of  absorbing  sphere  in  milk.  Axes  are  scaled  in  millimeters. 


Back-projected  image  presented  on  the  left  side  of  Figure  8  is  noisy.  Filtered  image  shown  on  the  right  side  of  the 
Fig.  8  demonstrates  the  noise  reduction.  The  contrast  of  filtered  image  (right)  was  increased  1.7  times  compared  with 
unfiltered  image  (left).  The  image  of  sphere  differs  from  the  correct  spherical  shape.  This  fact  has  a  simple  explanation.  As 

the  angle  of  acceptance  of  sphere  by  the  transducer  array  was  about  120°  only  one-third  of  data  needed  for  complete  image 
reconstruction  were  collected  and  employed  in  image  reconstruction. 


Figure  9.  Sphere  of  8-mm  diameter  in  milk  in  near  array  zone.  Milk  surface  is  marked  by  a  straight  line. 

An  increase  of  the  acceptance  angle  for  the  spherical  object  resulted  in  a  better  reproduction  of  its  shape  in  image,  as 
clearly  depicted  in  Fig.  9  showing  an  image  of  a  sphere  in  the  near  zone.  The  angle  of  acceptance  in  this  case  is  about  200° 
and  as  a  consequence  the  shape  of  image  highly  correlates  with  real  shape  of  the  sphere.  Fig.  9  (right)  presents  the  twice- 
magnified  fragment  of  the  sphere  image.  An  arc-shaped  artifact  could  be  seen  in  Fig.  9  and  resulted  from  the  opto-acoustic 
signal  generated  at  the  milk  surface.  This  surface  signal  was  not  filtered  completely  for  only  two  of  32  transducers,  but  it  was 
enough  to  yield  such  an  artifact.  Fig.  9  demonstrates  the  importance  of  optimal  signal  filtering  for  rejection  of  artifacts.  The 
filtration  procedure  demonstrated  in  Fig. 8  could  yield  complete  elimination  of  all  artifacts  due  to  surface  signals. 


3.2.  Spheres  in  gelatin  phantom 

A  gelatin  phantom  with  two  absorbing  spheres  was  made.  A  cylindrical  shape  of  the  phantom  with  120-mm 
diameter  provided  a  perfect  acoustics  contact  with  transducer  array.  The  effective  light  attenuation  coefficient  of  1.2  cm 
was  similar  to  typical  value  for  normal  breast  tissue  in  vivo  at  the  wavelength  of  1064  nm.  A  tumor  phantom  was  made  of  the 
same  gel  and  colored  with  the  bovine  hemoglobin.  The  “tumors”  had  spherical  shape  and  absorption  coefficient  of  1.0  cm*1. 
Diameter  of  the  spheres  was  approximately  8  mm.  The  relative  position  of  spheres  in  the  phantom  is  shown  schematically  in 
Fig.  10  (right  diagram).  The  depths  of  spheres  were  27  and  34  mm  respectively,  and  the  distance  between  their  centers  was 
about  1 1  mm.  The  procedure  of  measurements  was  as  follows. 


Laser  radiation 


Figure  10.  Image  of  two  spheres  in  gel  phantom. 

The  gelatin  phantom  was  placed  on  the  transducer  array  surface  for  the  opto-acoustic  scanning.  Absorbing  spheres 
were  located  above  the  transducer  array  beyond  its  focal  zone.  The  transducer  array  surface  was  preliminary  moistened  in 
order  to  provide  good  acoustic  coupling  between  phantom  and  the  transducer.  Optical  fiber  was  installed  at  the  10-mm 
height  from  the  phantom  surface.  The  fiber  was  moved  along  the  acoustic  array  axis  with  6-mm  intervals.  The  length  of  the 
opto-acoustic  scan  of  30-mm  was  sufficient  to  achieve  clear  resolution  of  images  of  each  sphere.  Full  set  of  32  acoustic 
signals  w  as  obtained  at  each  fiber  position.  As  a  result,  complete  table  of  data  corresponding  to  one  scan  was  acquired. 
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Figure  11.  OA  signal  of  two  absorbing  sphere  in  phantom  detected  by  one  of  the  central  transducer  in  array. 
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An  image  of  two  spheres  is  presented  in  Fig.  10  (left).  Procedure  of  image  filtering  was  employed  to  provide  the 
best  contrast  of  the  image  (-300%).  Two  spheres  are  clearly  resolved  and  depicted  as  two  separate  objects.  The  relative 
position  of  spheres  and  their  depths  are  in  good  agreement  with  their  real  location.  A  part  of  the  arc  phantom  surface  in  the 
poits  of  surface  illuminations  is  also  presented  in  Fig.  10.  The  image  of  the  surface  was  reconstructed  separately  using  the 
surface  marks  and  then  it  was  superimposed  on  the  sphere  image.  Partial  profile  of  a  signal  detected  by  one  of  the  central 
transducers  in  array  is  presented  in  Fig.ll.  Two  N-shaped  signals  that  were  excited  in  spheres  can  be  clearly  seen. 
Amplitude  of  the  signal  in  sphere  located  close  to  the  phantom  surface  (sphere  #1)  is  higher  then  in  sphere  #2.  Therefore,  the 
brigthness  and  the  image  contrast  is  higher  for  the  sphere  #1. 

3.3.  Breast  mastectomy  specimen 

Breast  specimens  were  obtained  after  radical  mastectomy  and  examined  within  the  first  ten  minutes  after  surgical 
excision.  It  provided  specimens  that  underwent  minimal  alterations  of  tissue  optical  properties,  except  some  blood  was 
drained  form  the  specimens.  The  mastectomy  specimens  were  measured  to  determine  tumor  geometry.  Location  of  tumors 
was  not  known  prior  to  opto-acoustic  imaging  procedure.  However,  tumors  were  sometimes  palpable  and  biopsy  incision 
was  visible  on  the  spared  segment  of  skin.  Opto-acoustic  image  of  breast  with  carcinoma  is  presented  in  Fig. 12. 
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Figure  12.  Opto-acoustic  image  of  a  malignant  tumor.  UTMB  Patient  #  4081 1 2P. 


Figure  13.  Ultrasound  image  of  tumor  UTMB  Patient  #  4081 12P. 
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Fifteen  fiber  positions  with  4-mm  interval  were  employed  for  data  collection  and  subsequent  image  reconstruction. 
The  scan  passed  directly  over  the  tumor.  The  tumor  core,  horizontal  tail  and  a  small  sprouts  may  be  clearly  seen  in  the  opto- 
acoustic  image  (confirmed  by  subsequent  pathology  examination).  The  dimension  of  a  tumor  core  can  be  estimated  as  10x7 
mm.  Tumor  was  detected  at  the  depth  of  1 1  mm  from  the  surface.  Skin  was  removed  from  most  part  of  the  specimen, 
therefore  the  irradiated  surface  is  depicted  as  a  rough  line.  Contrast  between  the  tumor  and  surrounding  normal  tissues  in  this 
image  exceeds  350%. 

The  ultrasound  image  of  the  same  tumor,  however  measured  in  vivo  is  presented  in  Fig.  13.  The  core  of  tumor  and 
its  long  tail  may  also  be  seen  in  this  image.  However,  small  sprouts  are  not  visible  to  limited  acoustic  contrast. 

4.  CONCLUSION 

The  second  generation  of  the  laser  optoacoustic  imaging  system  LOIS-02  for  breast  cancer  detection  with  32- 
element  transducer  array  was  developed  and  tested  in  phantoms  and  mastectomy  specimens.  Sensitivity  of  LOIS -02  permits 
detection  of  2-mm  tumors  at  the  depth  of  6  cm.  In  depth  resolution  equals  0.4  mm,  lateral  resolution  is  1-2  mm  depending  on 
position  of  tumor  relative  to  the  transducer  array. 

LOIS-02  is  the  basis  for  a  new  system,  which  will  be  used  for  in  vivo  studies.  In  vivo  optoacoustic  imaging  protocol 
is  currently  approved  at  UTMB.  Two  in  vivo  systems  are  being  fabricated:  compact  2D  real-time  system  and  3D  high- 
resolution  system  with  60-sec  data  acquisition  time. 
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ABSTRACT 

First  clinical  prototype  laser  optoacoustic  imaging  system  (LOIS)  for  breast  cancer  detection  was 
designed  and  fabricated  using  a  compact  Nd:YAG  laser,  fiberoptic  light  delivery  system,  a  linear  array  of  12 
wide-band  acoustic  transducers,  and  a  data  acquisition  card  operated  by  computer  with  original  signal  processing 
and  image  reconstruction  code.  Initially  images  of  small  absorbing  spheres  were  recorded  in  the  milk  with  optical 
properties  resembling  those  of  the  breast  at  the  wavelength  of  1064-nm.  The  system  was  optimized  for  contrast, 
sensitivity  and  axial  (in-depth)  resolution.  The  small  number  of  acoustic  transducers  (12),  which  in  turn  was 
determined  by  the  system  cost  and  the  time  of  image  acquisition  limited  the  lateral  resolution  of  the  images. 
Clinical  ex-vivo  studies  on  radical  mastectomy  specimens  were  performed  and  compared  with  x-ray  radiography, 
MRI  and  ultrasound  imaging.  The  results  of  our  pilot  clinical  studies  showed  pronounced  opto-acoustic  contrast 
of  -300%  between  breast  tumors  and  normal  breast  tissues.  This  contrast  substantially  exceeds  any  other 
endogenous  tissue  contrast  currently  utilized  in  clinical  ultrasonography,  MRI  and  x-ray  mammography.  Based 
on  literature  data  and  our  gross  observations  of  tumor  cross-sections  we  hypothesize  that  the  opto-acoustic 
contrast  results  primarily  from  increased  optical  absorption  in  the  dense  microvascularity  of  the  tumors.  In 
patients  receiving  radiotherapy,  tumors  were  found  to  contain  enhances  concentration  of  dense  highly  scattering 
fiberotic  tissue. 


1.  INTRODUCTION 

Breast  cancer  has  the  highest  incidence  (30%)  of  all  female  cancer  cases  in  the  USA.  About  10%  of 
American  women  develop  breast  cancer  during  their  lifetime.  It  is  estimated  that  175,000  new  cases  will  be 
diagnosed  and  43,300  deaths  will  occur  as  a  result  of  breast  cancer  this  year  [1].  Because  of  the  high  incidence, 
even  small  improvements  in  breast  cancer  diagnosis  and  treatment  may  save  tens  of  thousand  lives  every  year. 
Statistics  indicates  that  detection  and  follow-up  in  the  early  stages  of  cancer  offer  a  much  greater  chance  of 
survival  than  detection  and  follow-up  at  later  stages.  This  is  particularly  true  for  breast  cancer. 

The  gold  standard,  X-ray  radiography  (mammography),  magnetic  resonance  imaging  (MRI),  and 
ultrasound  are  the  three  imaging  techniques  that  are  currently  available  for  breast  cancer  detection.  All  these 
techniques  have  certain  advantages  for  detection  and  localization  of  a  specific  tissue  abnormality.  However,  all  of 
them  are  not  free  of  limitations  when  applied  for  the  breast  cancer  imaging.  The  major  limitation  of  X-ray 
mammography  is  its  incapability  to  detect  lesions  in  radiologically  dense  breast.  The  minimal  size  of  tumors 
detectable  by  mammography  is  about  5-10  mm  if  a  tumor  contains  no  calcification  [2].  X-ray  mammography 
utilizes  harmful  ionizing  radiation. 

The  main  limitation  of  the  ultrasound  imaging  is  a  low  contrast  [3].  On  one  hand,  the  detection  of  tumors 
in  many  cases  is  difficult  due  to  the  low  contrast  in  acoustic  properties  between  the  tumors  and  normal  tissues. 
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On  the  other  hand,  heterogeneous  structure  of  the  breast  yields  complex  ultrasound  images  with  numerous 
artifacts  resulting  from  multiple  reflections  of  the  ultrasonic  waves  from  tissue  boundaries.  The  minimal  tumor 
size  of  detection  by  current  ultrasonic  imaging  systems  at  a  depth  of  5-6  cm  inside  biological  tissues  is  limited  to 
about  10  mm. 

The  major  merit  of  MRI  is  high-resolution  [4].  This  technique  is  capable  of  characterizing  the  tissue 
hydrogen  content  and  water  environment.  However,  the  sensitivity  of  MRI  in  many  instances  is  not  very 
impressive,  and  this  technology  is  also  very  expensive  [4].  The  gadolinium  contrast  agent  is  used  in  MRI  in 
combination  with  dynamic  perfusion  imaging  to  visualize  tissue  vascularity,  which  helps  tremendously  to 
improve  cancer  detectability  [5, 6]. 

Various  types  of  optical  tomography  have  been  proposed  for  cancer  diagnostics  and  intensively 
developed  during  the  last  decade.  Optical  tomography  is  based  on  differences  in  optical  properties  between  tumor 
and  normal  tissues.  The  major  difference  in  absorption  coefficient  results  from  higher  blood  content  in  tumors 
compared  with  normal  tissues  [7]  which  in  turn  is  associated  with  increased  vascularization  (angiogenesis)  in 
rapidly  growing  tumors  [8,  9].  Since  hemoglobin  of  blood  is  the  major  chromophore  in  biological  tissues  in  the 
visible  and  near  infrared  spectral  range  [10,  11]  malignant  tumors  have  higher  absorption  coefficient  than  normal 
tissues.  Recent  clinical  studies  that  employed  optical  imaging  techniques  demonstrated  existence  of  optical 
contrast  of  200-300%  between  normal  and  cancerous  tissues  in  the  breast  [12,  13].  However,  strong  light 
scattering  and  attenuation  in  biological  tissues  substantially  limit  resolution  and  depth  of  monitoring  for  optical 
tomography. 

Several  years  ago  a  pulsed  laser  optoacoustic  tomography,  combining  advantages  of  optical  contrast  and 
sensitive,  high-resolution  ultrasonic  detection  was  proposed  for  breast  cancer  imaging  [14-17].  The  optoacoustic 
tomography  in  forward  mode  (OAT-F)  utilizes  acoustic  signals  induced  by  laser  pulses  in  tumors  and  transmitted 
forward  along  the  axis  of  laser  irradiation  to  the  acoustic  transducer  [18-20].  Application  of  a  wide-band  acoustic 
detection  instead  of  detection  of  photons  in  OAT-F  helps  to  overcome  two  problems  associated  with  strong  light 
scattering  in  biological  tissues  and  improve  depth  of  monitoring  [20]  and  spatial  resolution  [21,  22].  The  profiles 
of  acoustic  waves  generated  under  irradiation  conditions  of  temporal  pressure  confinement  in  the  volume  of 
tumors  resemble  the  profile  of  absorbed  laser  energy  [23].  Tumors  with  dimensions  of  1-10-mm  irradiated  with 
laser  pulses  represent  themselves  as  sources  of  acoustic  waves  with  ultrasonic  frequencies  of  -1  MHz  to  -100 
kHz.  Such  ultrasonic  waves  can  propagate  in  biological  tissues  with  insignificant  attenuation  [24-26],  However, 
acoustic  diffraction  in  the  ultrasonic  frequency  range  <  1  MHz  is  pronounced  and  must  be  compensated  by 
integrating  the  measured  signal  over  the  entire  time-course  of  detected  signal  [27,  28]. 

The  objectives  of  the  present  study  were:  (1)  to  construct,  assemble  and  test  in  vitro  the  first  clinical 
prototype  system  for  breast  imaging,  (2)  to  demonstrate  detection  and  localization  of  breast  tumors  in  mastectomy 
specimens  previously  visualized  in  vivo  with  conventional  imaging  modalities,  (3)  to  compare  contrast  between 
tumors  and  normal  tissues  for  four  techniques  (OAT,  X-ray  mammography,  MRI  and  ultrasound  imaging). 


2.  MATERIALS  AND  METHODS 


2.1.  Linear  array  transducer 

A  specially  designed  linear  array  of  PVDF  acoustic  transducers  fabricated  by  LaserSonix  Technologies, 
Inc.  Houston,  Texas,  was  employed  in  clinical  experiments.  A  photograph  of  the  transducer  array  is  shown  in 
Figure  1.  The  array  had  twelve  rectangular  piezo-elements  of  12.5x4.3  in  dimensions  and  1-mm  distance  between 
elements.  The  ultrasonic  frequency  band  determined  by  the  resonance  frequency  of  element  was  50  kHz  to  3 
MHz  at  the  1/e  level  of  maximal  sensitivity  of  30  //V/Pa.  The  axial  (in-depth)  resolution  of  images  obtained  with 
the  opto-acoustic  tomography  was  1-mm  as  determined  by  the  ultrasound  band  of  transducer  detection  [29].  The 
lateral  resolution  of  the  acoustic  transducer  array  was  defined  by  the  total  length  of  the  array  and  the  dimensions 
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of  each  individual  transducer  [29].  It  was  also  dependent  on  the  distance  between  the  tumor  and  the  transducer. 
An  average  lateral  resolution  in  our  clinical  experiments  was  approximately  equal  2a,  where  a  is  the  tumor  radius. 


Figure.  1.  Photograph  of  the  linear  array  of  12  transducers  housed  with  12  charge  preamplifiers. 

When  soft  tissue  with  light  absorbing  tumors  is  irradiated  with  a  10-nanosecond  laser  pulse,  it  is  heated 
faster  than  any  thermoelastic  expansion  may  occur.  Spherical  acoustic  waves  that  propagate  from  the  laser-heated 
tumors  have  bipolar  N-shaped  profiles.  The  duration  of  acoustic  profiles  excited  in  a  spherical  absorber  of  radius 
a0  equals  T  =  2Mq  /  Cs ,  where  cs  is  the  speed  of  sound.  The  smaller  the  tumor  size,  the  shorter  the  duration  of 
the  acoustic  N-pulse.  Wide  band  ultrasonic  frequency  detection  is  required  for  the  precise  transduction  of  the 
acoustic  profiles  into  electric  pulses.  A  two-dimensional  image  reconstruction  algorithm  employs  electric  pulses 
recorded  by  an  electronic  signal  acquisition  board  as  input  data. 

2.2.  Two  Dimensional  Opto-Acoustic  Imaging  System 

A  schematic  diagram  of  the  first  clinical  prototype  LOIS  is  shown  in  Fig.  2.  An  Nd:YAG  laser  operating 
at  the  wavelength  of  1064  nm  was  used  as  a  source  of  near-infra-red  pulses  of  10-ns  in  duration.  An  optical  fiber 
was  used  for  the  laser  pulse  delivery  to  the  tissue  surface.  The  fiber-optic  system  was  attached  to  a  three- 
dimensional  positioning  stage  driven  by  a  step  motor.  An  electronic  data  acquisition  system  was  synchronized 
with  the  linear  translation  along  the  chosen  direction.  Tissue  samples  were  placed  on  the  surface  of  a  50-mm 
thick  acoustic  buffer  filled  by  degassed  water.  The  buffer  thickness  was  sufficient  to  bring  the  rear  tissue  surface 
in  the  far  detection  zone.  Thereby  all  twelve  transducers  were  involved  in  the  opto-acoustic  detection,  which 
provided  a  better  lateral  resolution.  Amplified  electrical  signals  from  acoustic  transducers  were  digitized  by  a 
PDA-12  data  acquisition  card  with  two  input  channels  (SignaTech,  CA).  A  multiplexor  sequentially  connected 
twelve  amplifier  outputs  with  the  two  input  channels  of  the  data  card.  A  photodiode  detecting  laser  pulses  leaking 
from  the  rear  mirror  was  applied  for  synchronization  of  the  optical  and  acoustic  systems.  Acquired  data  were 
stored  in  a  computer  memory  for  further  processing.  The  data  processing  procedure  and  the  image  reconstruction 
algorithm  were  described  in  our  technical  publication  [29]. 

LOIS  sensitivity  is  a  function  of  the  system  detection  parameters  and  tissue  properties.  The  smallest  size 
of  detectable  absorbing  tumor  located  at  a  certain  depth  within  the  tissue  or  the  largest  depth  of  detection  for  a 
tumor  particle  with  certain  dimensions  can  characterize  the  sensitivity  of  the  opto-acoustic  imaging  system.  The 
light  absorption  coefficient  for  breast  tumors  depends  on  the  tumor  blood  content,  which  may  vary  in  a  wide 
range  around  ja^ 0.75  cm1  at  the  wavelength  of  1064-nm.  The  minimum  detectable  pressure  amplitude  is  limited 
by  the  thermoelectric  noise  of  the  piezoelement.  In  our  acoustic  transducers  the  effective  noise  pressure  may  be 
estimated  at  =0.5  Pa. 

In  order  to  increase  the  signal-to-noise  ratio,  an  increase  in  the  incident  laser  fluence  may  be  used.  However,  in 
clinical  applications  the  incident  laser  fluence  is  limited  by  a  safe  level  of  0. 1  J/cm2  permitted  by  the  American 


national  laser  safety  standards  [30].  Earlier,  we  demonstrated  experimentally  that  a  small  sphere  with  low  blood 
content  (/*a= 0.5  cm*1)  and  1-mm  radius  can  be  detected  at  the  depth  of  60  mm  when  a  safe  level  of  pulsed  laser 
irradiation  is  employed  [19]. 


Figure  2.  Schematic  diagram  of  LOIS. 

In  the  clinic,  the  breast  samples  were  examined  as  fast  as  possible  to  avoid  irreversible  chemical 
transformations  in  excised  tissue.  The  factor  of  time  is  critical  for  the  use  of  the  genetic  biochemical  assay  in 
pathology  studies.  The  system  was  designed  to  perform  opto-acoustic  imaging  procedure  in  less  than  10  minutes. 

2.3.  Conventional  imaging  modalities 

X-ray  images  of  the  breast  phantoms  were  obtained  with  the  use  of  clinical  mammography  machine 
(Lorad  M  IV,  Lorad  Corp.,  Danbury,  CT).  The  best  contrast  of  x-ray  images  was  achieved  with  application  of 
molybdenum  anode  and  molybdenum  filter  at  the  voltage  of  30  kV.  Ultrasonography  was  performed  with  a 
modern  sophisticated  ultrasound  machine  Sonoline  Elegra  (Siemens,  New  York,  NY)  with  various  contrast 
enhancement  features.  MRI  mammography  system  with  0.5  Tesla  magnet  (Caprius  Corp.,  Wilmington,  MA) 
employed  dynamic  perfusion  imaging  with  gadodiomide  contrast  agent  (Omniscan,  Nycomed-Amersham, 
Wayne,  PA). 

2.4.  Phantoms 

Breast  phantoms  were  made  of  10%-gelatin  and  had  dimensions  of  100  x  97  x  57  mm.  Water  absorption 
( jaa  =  0. 15  cm'1  at  1064  nm  [20])  provided  the  absorption  coefficient  of  0. 13  cm4  in  the  gel  phantoms.  The  typical 
effective  attenuation  coefficient,  piefU  for  normal  breast  tissue  in  vivo  is  about  1.0  to  1.2  cm'1  at  the  wavelength  of 
1064  nm,  while  the  optical  attenuation  of  the  skin  is  1.3- 1.4  cm4.  Polystyrene  spheres  were  used,  to  provide 
effective  scattering  coefficient  of  4  cm4.  Absorption  by  water  and  scattering  by  polystyrene  microspheres  yielded 
an  effective  attenuation  of  1.3  cm4  for  our  phantom. 

A  tumor  phantom  was  made  of  the  same  optically  turbid  gel  and  colored  with  the  bovine  hemoglobin. 
The  “tumor”  had  spherical  shape  and  absorption  coefficient  of  0.75  cm4.  Diameter  of  the  sphere  was 
approximately  8  mm.  The  sphere  was  placed  on  a  thin  nylon  white  thread  in  order  to  provide  its  fixed  position 
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inside  the  gel  phantom.  Liquid  gel  at  the  temperatures  very  close  to  the  solidification  temperature  was  poured 
inside  a  Plexiglas  cube.  After  solidification  in  a  refrigerator  the  phantom  was  taken  out  from  the  Plexiglas  cube 
and  examined  with  LOIS.  The  breast  phantoms  were  made  acoustically  homogeneous  and  radiologically  dense, 
while  the  ratio  of  absorption  coefficients  of  the  spheres  and  surrounding  gel  was  equal  6,  a  value  anticipated 
between  breast  tissues  and  advanced  tumors  in  vivo. 

2.5.  Mastectomy  specimen  and  the  imaging  procedure 

Breast  samples  were  obtained  after  radical  mastectomy  and  examined  within  first  ten  to  twenty  minutes 
after  excision.  It  provided  us  with  specimen  that  underwent  minimal  alterations  of  tissue  optical  properties, 
except  the  blood  was  noticeably  drained  form  the  specimen.  The  mastectomy  specimens  were  measured  to 
determine  their  dimensions.  Location  of  tumors  was  not  known  prior  to  opto-acoustic  imaging  procedure. 
However,  often  tumors  were  palpable  and  a  biopsy  incision  was  visible  on  the  spared  segment  of  the  skin. 


Figure  3.  Photograph  of  mastectomy  specimen  placed  on  the  transducer  array  for  imaging  with  LOIS. 

Breast  samples  were  placed  on  the  transducer  buffer  surface  for  the  optoacoustic  scanning.  The  buffer 
surface  was  preliminary  wetted  to  provide  a  good  acoustic  contact  between  breast  tissues  and  the  transducer. 
Optica]  fiber  was  installed  at  the  3-mm  height  from  the  sample  surface.  The  fiber  was  moved  along  the  acoustic 
array  axis  with  6-mm  intervals.  The  laser  beam  diameter  on  tissue  surface  was  about  8-mm.  The  repetition  rate 
was  10  Hz  and  the  pulsed  energy  was  approximately  48  mJ.  The  length  of  the  opto-acoustic  scan  of  72-mm  was 
equal  to  the  length  of  the  acoustic  transducer  array.  At  each  fiber  position  a  set  of  acoustic  signals  were  obtained. 
As  a  result  the  complete  table  of  data  corresponded  to  one  scan  was  acquired.  Imaging  procedure  (data 
acquisition)  was  performed  in  vertical  direction  (along  the  body  axis)  two-three  times  within  portions  of  the  breast 
suspicious  to  have  a  tumor. 

Gross  photographs  of  the  specimens  were  taken  before  and  during  pathological  examination.  The  first 
stage  of  the  pathological  examination  (performed  immediately  after  the  opto-acoustic  tomography)  included 
sectioning  of  the  breast  specimen  in  vertical  direction  producing  0.5-cm  to  1-cm  thick  slices.  Opto-acoustic 
images  were  reconstructed  based  on  the  measured  pressure  profiles.  The  images  were  compared  with  x-ray 
mammograms,  ultrasonograms  and  MRI  images,  when  available,  and  with  gross  and  histological  pathology 
findings  obtained  during  routine  pathologic  examination. 


3.  RESULTS 

A  typical  image  of  a  sphere  8-mm  in  diameter,  embedded  in  gel  phantom  is  shown  in  Figure  3.  The 
center  of  sphere  was  located  at  the  50-mm  depth  from  the  irradiated  surface,  and  its  center  was  shifted  8-mm  to 
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the  left  from  the  phantom’s  center.  This  model  experiment  was  performed  to  test  the  sensitivity  and  resolution  of 
our  optoacoustic  imaging  system. 


Figure  3.  A  8-mm  diameter  gel  sphere  located  at  a  depth  of  5-cm  in  milky  phantom  resembling  optical  properties 
of  the  breast.  Axes  are  scaled  in  millimeters.  The  optical  contrast  of  the  sphere  above  the  background 
medium  equals  5. 

The  top  of  Fig.  3  corresponds  to  the  location  of  optical  fibers  (Y=0  mm),  while  bottom  of  this  figure 
corresponds  to  the  position  of  transducer  array  (Y=l  10  mm).  The  center  of  the  transducer  array  was  at  the  zero  of 
X-axis.  The  opto-acoustic  image  in  Fig.  3  represents  the  distribution  of  a  product  of  thermo-acoustic  efficiency, 
optical  absorption  coefficient  and  the  absorbed  laser  energy.  The  brightest  parts  of  the  image  correspond  to  the 
maximum  values  of  this  product.  The  sphere  location  and  its  vertical  dimension  accurately  correspond  to  their 
real  values.  The  lateral  size  of  image  is  approximately  3  times  larger  than  the  sphere  diameter.  It  demonstrates 
the  limit  of  lateral  resolution  of  the  system.  Five  positions  of  illuminating  Fiber  with  6-mm  increment  were 
employed  to  collect  data  for  this  image  reconstruction. 


Figure  4.  Opto-acoustic  image  of  a  large  tumor.  Scan  #3.  Patient  052291P. 


Typical  results  of  the  clinical  opto-acoustic  imaging  are  presented  in  Figures  4  -  9. 


The  laser,  the  acoustic  transducer  array  and  the  data  acquisition  system  were  placed  on  a  cart.  The  cart 
provided  means  for  fast  transportation  of  the  system  to  the  frozen  section  room  in  Pathology  Department  where 
all  measurements  were  performed. 

The  image  of  large  well  palpable  tumor  is  shown  in  Fig.  4.  Total  of  four  opto- acoustic  scans  were  made. 
The  data  for  the  image  presented  in  Fig.  4  were  obtained  when  the  scan  passed  directly  over  the  tumor  (scan  #3). 
Tumor  position  and  direction  of  scans  are  shown  in  the  diagram  (on  the  right).  Axes  in  Fig.  4  are  scaled  in 
millimeters.  The  top  of  the  figure  (Y  =  0  mm)  corresponds  to  the  skin  surface,  while  the  transducer  array  was 
located  on  the  bottom  (Y  =  80  mm).  Two  bright  domains  in  Fig.  4  depict  a  two-lump  tumor  (a  smaller  lump  was 
located  underneath  the  larger  lump.  The  full  vertical  size  of  this  tumor  could  be  determined  as  16-18  mm.  The 
depth  of  the  tumor  was  found  to  be  10-12  mm  from  the  skin  surface.  Spherical  shadow  arcs  in  the  tumor  image 
appeared  due  to  the  secondary  lobes  in  the  sensitivity  diagram  of  the  ultrasonic  reception  of  the  array  transducer 
[29].  This  type  reception  diagram  represents  the  major  limitation  for  the  lateral  resolution  in  our  imaging  system. 

The  temporal  integral  of  the  initial  signal,  detected  by  one  of  the  transducers  in  array  is  shown  in  Fig.  5. 
Initial  signal  consists  of  the  three  components:  short  N-shaped  signals  generated  by  the  small  absorbing  tumors,  a 
lower-frequency  smooth-rise  exponential  slope  produced  by  the  optical  attenuation  inside  the  tissue,  and 
transducer  reverberations  detected  after  the  end  of  useful  part  of  the  signal. 


Time,  ps 

Figure  5.  Integrated  opto-acoustic  signal  profiles  measured  in  the  mastectomy  specimen  with  tumor. 

It  is  necessary  to  eliminate  completely  the  low  frequency  signal  component  associated  with  homogeneous 
light  attenuation  inside  the  tissue,  because  it  can  significantly  decrease  the  contrast  in  the  process  of  image 
reconstruction.  This  was  achieved  in  LOIS  with  a  high-pass  digital  filtering  using  RC  cells.  Variation  of  the  RC 
parameters  and  the  number  of  cells  permits  a  fine  adjustment  of  the  filter  bandwidth  and  sharpness  of  this  band 
profile. 


High  ultrasound  frequencies  that  fall  near  the  resonance  frequency,  /max,  of  piezo-elements  can  induce 
significant  reverberation  in  the  transducer.  Our  transducers  were  designed  to  damp  resonance  frequencies  and 
widen  the  range  of  detectable  frequencies.  Nevertheless,  a  strong  reverberation  takes  place  after  the  arrival  of  the 
signal  from  the  surface  of  tissue.  These  reverberations  do  not  affect  the  useful  part  of  the  opto-acoustic  profile. 
Besides,  it  is  convenient  to  employ  the  sharp  gradient  in  the  signal  for  automatic  determination  of  the  tissue 
surface.  The  position  of  maximum  in  the  signal’s  derivative  was  considered  as  the  position  of  the  irradiated  tissue 
surface.  All  data  samples  were  cut  off  the  signal  after  arrival  of  the  surface  signal.  Nevertheless,  the  signal 
amplitude  near  the  surface  position  remains  high  and  could  give  rise  to  a  bright  area  on  images.  To  avoid  the 
surface  signal  and  to  enhance  the  tumor  image  contrast  the  surface  signal  was  also  rejected. 
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The  image  of  the  same  tumor  but  for  another  scan  direction  (scan#4)  is  shown  in  Fig.  6.  This  scan  was 
made  aside  the  tumor  center.  As  a  result,  the  depth  of  tumor  image  is  slightly  increased  (14  -15  mm)  and  its 
vertical  size  is  decreased  (13-14  mm)  compared  with  Fig.  4. 


Figure  6.  Image  of  tumor.  Scan  #4.  ITMH  Patient  #052291P. 

The  main  tumor  lump  (the  brightest  part  on  the  image)  is  shifted  towards  the  edge  of  the  transducer  array. 
At  the  depth  of  40  mm  a  layer  of  blood  on  the  surface  of  the  transducer  buffer  is  visualized.  Blood  was  mixed 
with  water  and  this  solution  was  accumulated  at  the  surface  of  the  film  covered  the  transducer  buffer.  The  light 
absorption  in  this  layer  resulted  in  additional  acoustic  signal  and  this  layer  was  imaged.  The  depth  of  this  layer  of 
diluted  blood  was  relatively  large  (40  mm  from  the  illuminating  surface)  and  the  blood  concentration  was  weak. 
Thus,  the  visualization  of  this  layer  demonstrated  high  sensitivity  of  LOIS.  As  it  was  mentioned  earlier  the  tumor 
depicted  in  Fig.  4,  and  Fig.  5  was  large  and  it  was  detected  and  visualized  as  well  by  other  imaging  modalities  (X- 
ray  mammography,  ultrasound,  and  MRI). 


Figure  7.  Opto-acoustic  signal  measured  in  normal  breast  tissue. 

Fig.  7  presents  acoustic  signals  exited  in  normal  breast  tissue.  The  acoustic  signal  (dotted  curve)  has  a 
smooth  shape  without  any  peaks  associated  with  additional  light  absorption.  This  form  can  be  approximated  by 
exponential  function  in  order  to  determine  the  coefficient  of  optical  attenuation  in  the  breast.  The  optical 
attenuation  of  the  normal  breast  tissues  at  the  wavelength  of  1064  nm  was  determined  from  Fig.  6  to  be  \ic»  —  1.2 
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cm  l.  When  the  same  tissue  was  illuminated  through  the  skin  layer  (solid  curve),  an  additional  light  scattering  and 
attenuation  in  skin  resulted  in  two-peak  profile  with  increased  amplitude  of  the  acoustic  profile. 

The  image  of  tumor  attached  to  the  surface  of  silicon  implant  is  presented  in  MRI  image  and  ultrasound 
image  depicted  in  Fig.  8.  The  presence  of  implant  complicated  all  imaging  procedures,  resulting  in  a  very  low 
contrast  between  the  tumor  and  the  implant.  X-ray  mammography  (not  shown)  did  not  reveal  any  tumors.  Both 
MRI  and  ultrasound  depicted  area  suspicious  to  have  a  tumor  (see  a  brighter  area  on  the  top  surface  of  the 
implant).  MRI  system  employed  dynamic  perfusion  imaging  with  gadodiomide  contrast  agent  that  resulted  in  a 
better  differentiation  between  tumor  and  the  implant.  However,  the  contrast  between  the  tumor  and  surrounding 
tissues  was  not  sufficient  to  make  pathological  conclusions.  A  biopsy  performed  following  imaging  demonstrated 
the  presence  of  malignant  cells  in  the  suspicious  area  on  implant  surface. 


Figure  8.  Differential  MRI  image  of  the  breast  with  tumor  (left,  1  minute  after  contrast  minus  pre-contrast),  and 
Ultrasound  image  of  the  same  breast,  employing  7.5  MHz  transducer. 


Figure  9.  Tumor  on  the  implant  surface.  UTMB  Patient  #166248P. 
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The  corresponding  opto-acoustic  image  is  shown  in  Figure  9.  The  initial  raw  data  for  this  image 
reconstruction  were  acquired  during  the  scan  shown  in  the  diagram  (to  the  right  from  the  image).  The  implant 
was  about  40  mm  in  thickness.  The  tumor  (bright  domain  at  the  depth  of  about  30-mm  from  the  skin  surface)  was 
detected  near  the  top  edge  of  implant.  It  is  important  to  note  that  the  tumor  location  on  the  opto-acoustic  image 
differs  from  that  on  the  MRI  image  to  the  factors  influencing  imaging  of  mastectomy  specimen,  i.e.  (I)  different 
shape  and  dimensions  of  the  breast  in  vivo  and  after  surgical  excision,  and  (2)  shifted  position  of  the  mastectomy 
specimen  between  the  fiber  and  the  transducer  array  compared  with  the  natural  position.  The  area  where  the 
enhanced  light  absorption  was  found  had  5-7  mm  thickness.  Some  blood  drained  to  the  transducer  buffer  surface 
(at  the  depth  of  approximately  75  mm)  and  mixed  with  water  resulted  in  a  visible  mark  on  the  image.  The  opto- 
acoustic  tissue  contrast  between  the  tumor  and  the  surrounding  tissue  and  the  implant  is  significant.  However,  an 
enhanced  contrast  could  result  from  distortion  of  breast  structure  and  blood  accumulation  after  the  biopsy. 

4.  DISCUSSION 

This  comparative  study  demonstrated  that  the  laser  opto-acoustic  imaging  could  provide  images  with 
substantially  better  contrast  than  x-ray  mammography  and  ultrasound  imaging  in  case  of  radiologically  dense, 
acoustically  homogenous,  but  optically  inhomogeneous  phantoms.  Contrast  and  sensitivity  are  the  two 
interdependent  parameters.  The  higher  the  tissue  contrast,  the  higher  the  sensitivity  of  the  opto-acoustic  system  in 
detecting  small  tumors.  Effective  noise  pressure  could  be  determined  from  the  following  expression: 

,  2p0c5/max  [4kT 

Pn= - (0 

833  1  Cq 

where  g33  is  the  piezoelectric  modulus,  fmax  is  the  maximum  detected  ultrasound  frequency,  C0  is  the  transducer 
capacity. 

Figure  10  shows  effective  noise  pressure  as  a  function  of  the  transducer  capacity  for  gradually  increasing 
width  of  detected  ultrasonic  frequencies. 


Piezoetement  Capacity,  pF 

Figure  8.  Effective  noise  pressure  of  the  acoustic  transducer. 

A  typical  effective  noise  pressure  is  about  10  /*bar=l  Pa  for  the  acoustic  transducers  with  capacity  of 
about  10  pF. 
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The  axial  (in-depth)  resolution  of  the  pulsed  opto-acoustic  system  operating  under  irradiation  conditions 
of  temporal  pressure  confinement  is  directly  defined  by  the  detection  bandwidth  of  acoustic  transducers.  Two 
coherently  emitting  spherical  acoustic  sources  with  radius  a0  will  be  detected  as  separate  objects  if  the  space 
between  their  centers  equals  or  greater  than  Az  =  2a0.  For  correct  detection  of  pulsed  opto-acoustic  waveform 
excited  in  a  tumor  with  characteristic  dimension  a0  it  is  necessary  to  detect  opto-acoustic  signals  within  a  wide 
ultrasound  frequency  band  equal  approximately  the  ratio  of  the  Cs /  a0.  Lateral  resolution,  <5,  provided  by  an  array 
of  transducers  can  be  estimated  as  a  product  of  the  angle  of  acceptance  (numerical  aperture)  of  the  acoustic 
transducer  and  the  total  length  of  the  transducer  array. 


5.  CONCLUSION 

The  study  performed  with  breast  phantoms  demonstrated  capabilities  of  the  laser  opto-acoustic 
tomography  in  detection  of  small  volumes  with  absorption  equal  to  the  absorption  of  10%  blood  at  the  depth  of  up 
to  6  centimeters  in  phantoms  resembling  optical  properties  of  breast  tissues.  Thus,  the  laser  optoacoustic 
tomography  has  a  potential  for  early  detection  of  (small)  breast  tumors. 

The  first  clinical  prototype  of  the  laser  opto-acoustic  imaging  system  was  developed  and  tested  on  radical 
mastectomy  specimen.  Significant  opto-acoustic  contrast  (-300%)  between  tumors  and  surrounding  tissues  was 
observed  and  exceeded  the  contrast  of  conventional  imaging  modalities.  LOIS  was  capable  of  detection  all  and 
every  tumor  in  mastectomy  specimens.  Because  modern  medical  diagnostic  techniques  have  limitations  in 
detection  of  small  tumors  (with  a  volume  of  <  4  mm3)  located  at  depths  of  5-6  centimeters,  the  laser  opto-acoustic 
imaging  may  occupy  a  significant  niche  in  diagnosis  of  breast  tumors,  especially  in  radiologically  dense  and 
acoustically  homogenous  breasts. 
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ABSTRACT 

The  major  limitation  in  sensitivity  of  the  optical  tomography  is  associated  with  strong  optical  attenuation  in  human  tissues. 
Opto-Acoustic  Tomography  overcomes  this  limitation  utilizing  detection  of  acoustic  waves  instead  of  detection  of 
transmitted  photons.  Exceptional  sensitivity  of  the  opto-acoustic  tomography  allows  early  detection  of  small  tumors  located 
deep  in  human  tissues,  such  as  breast.  This  paper  demonstrates  that  an  optimally  designed  opto-acoustic  imaging  system  can 
detect  early  I -mm  tumors  with  minimal  blood  content  of  only  7%  at  the  depth  of  up  to  7-cm  within  the  breast  attenuating 
laser  irradiation  3.3  times  per  each  1-cm  of  its  depth.  A  theoretical  consideration  of  the  ultimate  sensitivity  of  piezo¬ 
detection  in  a  wide  ultrasonic  frequency  band  is  developed.  The  detection  sensitivity  is  presented  as  a  function  of  the 
ultrasonic  frequency,  tumor  dimensions  and  optical  absorption  coefficient.  Comparative  analysis  of  piezo  and  optical 
interferometric  detection  of  opto-acoustic  transients  is  presented.  The  theoretical  models  of  piezodetection  were  developed 
for  the  open-circuit  and  short-circuit  schemes  of  operation.  The  ultimate  sensitivity  limited  by  thermal  noise  of  electric 
capacitor  of  the  piezo-element  was  estimated.  It  was  shown  that  the  limit  of  detection  depends  on  the  frequency  band,  the 
electric  capacity  of  the  transducer  and  the  speed  of  sound  in  the  piezo-element.  Comparative  analysis  of  various  piezo¬ 
materials  was  made  from  the  point  of  view  of  their  utility  for  sensitive  opto-acoustic  detection. 

Keywords:  piezoelectric  transducer,  optical  vibrometer,  opto-acoustic  imaging 


1.  INTRODUCTION 

Advantages  of  the  opto-acoustic  imaging  for  characterization  of  optically  heterogeneous  tissue  structures  with  high 
spatial  resolution  were  demonstrated  over  the  past  several  years  [1-16].  One  of  the  major  merits-of  the  opto-acoustic  imaging 
is  high  sensitivity  in  detecting  variations  of  laser-induced  temperature  or  corresponding  pressure.  Unique  characteristics  of 
optoacoustics  found  utility  in  various  medical  applications,  including  detection  of  early  cancer  [7-1 1],  monitoring  therapeutic 
effects  of  lasers  in  tissues  [12,13],  measuring  tissue  properties  and  molecular  content  [1,14-16].  Optimal  design  and 
development  of  optoacoustic  imaging  and  monitoring  systems  requires,  first  of  all,  comprehensive  analysis  of  capabilities 
and  limitations  in  detection  of  wide-band  ultrasonic  transients. 

There  are  two  general  modalities  that  may  be  used  for  conversion  of  wide-band  acoustic  waves  (mechanical  displacement) 
into  electronic  signal:  piezoelectric  transducers  [5,17-19]  and  optical  interferometers  and  deflectors  [20-22].  Noncontact 
nature  of  laser  generation  of  ultrasound  makes  it  attractive  to  develop  fully  noncontact  opto-acoustic  tomography  and  laser 
ultrasound  evaluation  systems  employing  optical  methods  of  detection  of  acoustic  transients.  However,  surface  roughness 
(especially  in  skin  and  other  uneven  media)  strongly  compromises  sensitivity  of  optical  detection.  The  use  any  contact 
means  for  surface  smoothing  removes  major  advantage  of  noncontact  optical  methods  relative  to  the  contact  piezoelectric 
detection.  Therefore,  the  main  attention  is  drawn  in  this  paper  to  the  piezoelectric  detection  of  acoustic  waves.  The  most 
sensitive  acoustic  detectors  utilize  resonance  of  physical  dimensions  and  ultrasound  wave  to  be  detected.  The  product  of 
sensitivity  and  the  bandwidth  is  constant  defined  by  properties  of  the  detector.  Therefore,  detection  within  a  narrow 
frequency  band  is  the  most  sensitive.  The  sensitivity  goes  down  simultaneously  with  widening  of  the  ultrasound  detection 
band.  In  order  to  account  for  the  fact  that  sensitivity  of  acoustic  transducer  depends  on  the  range  and  the  bandwidth  of 
ultrasound  detection,  one  should  consider  the  bandwidth-specific  sensitivity  of  the  ultrasonic  detector. 

The  aim  of  this  paper  is  to  describe  ultimate  sensitivity  of  detection  of  acoustic  waves  in  the  entire  ultrasonic 
frequency  band  from  20kHz  to  200  MHz  and  to  compare  piezoelectric  transducers  versus  optical  interferometry  as  means  for 
detecting  transient  ultrasonic  signals  induced  in  tissue  with  laser  pulses.  It  is  logical  to  divide  all  applications  of  the  opto¬ 
acoustic  imaging  in  three  groups  according  to  the  required  bandwidth  of  ultrasound  detection:  low,  intermediate  and  high,  (1) 
20  kHz  to  2  MHz,  (2)  1  MHz  to  20  MHz,  and  (3)  1  MHz  to  200  MHz. 
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The  first  group  includes  imaging  and  monitoring  of  tumors,  blood  vessels  and  other  tissues  located  at  the  depth  of  up 
to  10  cm.  Substantial  depth  of  imaging  prohibits  the  use  of  high  ultrasound  frequencies  due  to  their  strong  attenuation  in 
tissues.  Therefore,  spatial  resolution  is  limited  to  sub-mm.  The  second  group  includes  applications  of  tissue  imaging  at  the 
depth  of  about  1-cm.  These  applications  require  a  more  detailed  information  of  tissue  structure  and  function,  and  therefore, 
employ  detection  of  higher  ultrasound  frequencies  in  order  to  achieve  resolution  of  several  tens  of  microns.  The  third  group 
includes  applications  of  the  opto-acoustic  tomography  in  high  resolution  imaging  of  skin  and  mucosal  surfaces  of  digestive 
tract,  and  in  cell  microscopy.  Each  group  of  applications  and  each  corresponding  range  of  detected  ultrasound  frequencies 
define  a  unique  design  of  the  opto-acoustic  detection  system. 


2.  ULTRASONIC  SPECTRUM  OF  THE  OPTO-ACOUSTIC  SIGNAL 


Opto-acoustic  imaging  of  biological  tissues  utilizes  absorption  of  optical  energy  by  heterogeneous  structures  of 
tissue  and  materials  and  detection  with  temporal  resolution  of  the  wide-band  ultrasound  transients  resulting  from 
nonstationary  thermal  expansion  [23].  The  profile  of  acoustic  wave  emitted  by  a  source  with  complex  shape  may  be 
represented  by  the  convolution  of  acoustic  waves  emitted  by  a  number  of  small  spherical  sources.  Therefore,  below  acoustic 
emission  is  described  for  a  small  optically  absorbing  sphere  placed  in  nonabsorbing  medium  with  thermomechanical 
properties  similar  to  those  of  the  sphere.  Let  us  consider  the  case  when  laser-induced  heat  is  evenly  distributed  over  the 
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spherical  volume  of  — - —  ,  where  Rp  is  the  radius  of  the  sphere.  Other  types  of  radiation  (such  as  microwave,  IR  or  x- 
ray)  and  other  shapes  of  absorbing  object  will  not  change  our  model  qualitatively. 

Acoustic  pulse  generated  upon  instant  heating  of  the  spherical  volume  can  be  expressed  as  in  [24]: 


p'(T  =  t  -  r / cq)  = 
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C0T,  |c0r|  ^  R0 , 
0,  |c0t|  >  R0\ 


(1) 


where  p  (r)  is  the  pressure  increse  in  ultrasonic  wave  propagating  from  the  center  of  the  sphere,  (3  is  the  thermal 
expansion  coefficient,  Cp  is  the  specific  heat-capacity,  Eafrs  is  the  optical  energy  absorbed  by  the  sphere,  r  is  the  distance 
from  the  point  of  observation  to  the  center  of  the  spherical  particle. 


Retarded  time 


Figure  1.  Typical  N-shaped  laser-induced  acoustic  profile  of  from  a  small  absorbing  sphere  with  radium  Rp. 

The  opto-acoustic  profile  has  an  N-shape  profile  with  traingle  compression  phase  followed  by  the  antisimmetric 
tension  phase  (see  Figure  1).  The  total  duration  of  the  opto-acoustic  signal  is  defined  by  the  time  of  sound  propagation  along 
the  diameter  of  the  sphere.  Since  we  consider  the  bandwidth  specific  sensitivity  of  the  opto-acoustic  system,  let  us  consider 
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the  spectrum  of  the  ultrasonic  signal.  The  ultrasonic  spectrum  of  the  opto-acoustic  signal  (1)  could  be  obtained  taking 
Fourier  transform  of  the  signal's  temporal  profile. 
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The  resulting  spectrum  has  the  following  form: 


PaM)'* 


(coRp  /  CQ )  cos  (0)Rp  /  CQ )  -  sinjcoRp  /  c0 ) 
(coRp/c0)2 


where  is  the  frequency  of  ultrasonic  vibrations.  The  modulus  of  this  spectrum  is  given  in  Figure  2  as  a  sequence  of 
oscillations  rapidly  decreasing  in  their  amplitude.  Maximum  of  the  ultrasonic  emission  spectrum  is  given  by  the  frequency, 
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max 


=  0.33— 

R„ 


Normalized  frequency  (2 n  f  RQ/  cQ) 

Figure  2.  Fourier  spectrum  of  acoustic  emission  from  a  small  sphere  as  a  function  of  ultrasonic  frequencies. 
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The  bandwidth  of  this  spectrum  is  usually  measured  at  the  level  of  0.707  or  -3dB  between  —  0.16—  and 

Rp 

Q o 

f upper  =  0*5 1-~".  On  the  other  hand,  one  has  to  take  into  account  the  presence  of  various  absorbing  objects  in  the  laser- 

Rp 

irradiated  volume  and  detect  all  ultrasonic  frequencies  corresponding  to  variety  of  sphere  dimensions  from  their  minimum  to 
maximum.  Thus,  wide  band  detectors  in  which  the  upper  limit  of  the  ultrasonic  detection  band  exceeds  the  lower  limit  more 
than  10  times  are  required  for  correct  opto-acoustic  detection  of  laser-induced  acoustic  sources.  Below  we  discuss  noise 
level,  signal-to-noise  ratio  (SNR)  and  sensitivity  of  this  type  wide-band  acoustic  trasnducers. 


3.  PIEZOELECTRIC  DETECTION  OF  WIDE-BAND  ULTRASONIC  TRANSIENTS 

There  are  two  operation  modes  for  wide  band  piezoelectric  transducers:  short  circuit  and  open  circuit  (idling)  [25]. 
The  electrical  circuits  of  the  two  modes  of  transducer  operation  are  presented  in  Figures  3a  and  3b. 

The  transducer  operating  in  the  short  circuit  has  substantial  thickness  that  is  larger  than  the  spatial  width  of  the 
detected  ultrasonic  transient  (Figure  3a).  The  transducer  thickness  limits  the  lower  limit  of  detected  ultrasonic  frequencies 
and  duration  of  the  detection  gate.  An  excessively  big  thickness  of  the  piezoelement  would  (1)  lead  to  a  more  prominent 


acoustic  diffraction  at  lower  ultrasonic  frequencies  in  the  detected  signal  and  (2)  reduce  the  electric  capacity  to  the  value 
below  the  electric  capacity  of  the  electronic  circuitry.  It  is  difficult  to  design  acoustic  transducers  operating  in  short-circuit 
mode  for  detection  of  opto-acoustic  profiles  longer  than  1  ps.  However,  these  transducers  would  be  most  optimal  for 
detection  of  submicrosecond  and  nanosecond  opto-acoustic  signals.  The  upper  limit  of  ultrasonic  frequency  is  defined  by  the 
discharge  time  of  the  transducer  capacity,  Cj  through  the  resistor,  R.  With  low  values  of  RCt  the  ultrasonic  detection  band 
could  reach  several  hundred  MHz.  The  optimal  value  of  the  resistor  R  equals  wave  resistance  of  the  electrical  cable  (50  to 
200  Ohm). 


(a) 


- 1 - *-► 

OPEN  (IDLER)  CIRCUIT  MODE 

(b) 


Figure  3.  Electric  circutry  for  the  short  (a)  and  idler  (b)  operation  modes  of  the  wide-band  acoustic  transducer. 

In  case  of  acoustic  transducer  operating  in  idler  circuit  mode  (Figure  3b)  it  is  necessary  that  the  thickness  of 
piezoelectric  element  be  smaller  than  the  acoustic  wavelength  detected  in  its  upper  ultrasonic  frequency  limit.  The  lower 
limit  of  detectable  ultrasonic  frequencies  is  defined  by  the  discharge  time  of  transducer  electric  capacity,  Ct,  through  the 
input  resistor  of  the  electronic -preamplifier,  R.  Therefore,  it  is  desirable  to  design  these  acoustic  transducers  in  one  and  the 
same  housing  with  charge  preamplifier.  In  order  to  design  acoustic  transducers  operating  in  open  (idler)  circuit  mode  with 
very  wide  ultrasonic  detection  band  (>100  MHz),  one  have  to  employ  very  thin  piezoelectric  elements  (thickness  of  just  a 
few  microns).  Technocal  difficulties  are  associated  with  flatness  of  the  piezoelement  and  necessity  for  effective  damping  of 
high  frequency  resonances.  Design  of  the  backing  layer  acoustically  matched  with  piezoelectric  material  for  effecting 
damping  of  transducer  resonances  is  one  of  the  most  important  technical  problems.  These  technical  difficulties,  however,  are 
compemsated  with  a  higher  sensitivity  of  transducers  operating  in  the  open  circuit  mode.  The  sensitivity  of  open  circuit 
transducers  is  about  50  times  that  of  short  circuit  transducers  due  to  longer  holding  of  electric  charge  at  the  piezoelement. 


3.1.  Thermal  noise  in  acoustic  transducers  operating  in  the  open  circuit  mode. 

Let  us  consider  a  thin  piezoelectric  layer  placed  in  medium  with  ideally  matching  acoustic  impedance  and  operated 
in  open  electric  circuit  (Fig. 3b).  The  upper  limit  of  detectable  ultrasonic  frequency  determinde  at  the  level  of-3dB  may  be 
found  from  the  following  expression  [26]: 


fupper  =  0.44-^  ,  (5) 

where  l  is  the  thickness  of  piezoelectric  element,  C/  is  the  velocity  of  sound  in  piezoelectric  material. 

Therefore,  minimal  electric  capacity  of  piezoelectric  element  is  defined  by  the  upper  limit  of  ultrasonic  detection 

band  as: 


C  >  2-3££qA  ~ 

—  J upper  • 


(6) 


where  A  is  the  area  of  piezoelectric  element,  £  and  £0  are  the  dielectric  constants  of  piezoelectric  material  and  vacuum. 

Formula  (6)  shows  that  the  choice  of  piezoelectric  material  for  wide-band  acoustic  transducers  is  affected  not  only 
by  the  dielectric  permeability,  as  in  case  of  resonance  acoustic  transducers,  but  also  by  the  velocity  of  sound  in  piezoelectric 
material  (in  addition  to  upper  ultrasonic  frequency,  fUpper)> 
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For  a  transducer  operating  in  the  open  circuit  mode,  the  electrical  conductivity  of  piezoelectric  element  substantailly 
exceeds  the  electric  conductivity  of  the  load  resistor,  i.e.  the  current  through  piezoelement  is  greater  than  the  drain-current 
through  the  resistor.  Therefore,  the  voltage  due  to  thermal  noise  of  the  piezoelectric  element  of  transducer  can  be  found  from 
thermal  energy  considerations  (Neiquist  formula)  as: 


where  k  is  the  Boltsman  constant,  TQ  is  the  temperature  of  piezoelement. 
The  voltage  of  piezoelectric  effect  is  given  by: 


(7) 


Us  -  833I  p' . 


(8) 


where  g  is  the  piezoelectric  pressure  modulus  of  the  transducer  material.  Sometimes,  it  is  more  convenient  to  use 

2 

piezoelectric  deformation  modulus,  =  <?33  PQ  for  description  of  the  piezoelectric  sensitivity. 

The  pressure  corresponding  to  the  thermal  noise  voltage  can  be  expressed  combining  formulas  (7)  and  (8): 


\9-2kTof upper  _  1 92kTo 

V  cl££0§23A  V  Vac  A 

The  expression  (9)  indicates  that  the  level  of  pressure  minimally  detectable  with  wide-band  acoustic  transducers 
operating  in  open  circuit  mode  increases  with  increasing  ultrasonic  detection  bandwidth.  Here  we  introduced  a  parameter  of 
piezoelectric  element  important  for  wide-band  ultrasound  detection  and  termed  as  «bandwidth-specific  piezoelectric 
efficiency)): 
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This  parameter  defines  optimal  choice  of  the  materizl  for  sensitive  wide-band  ultrasonic  detection.  Formula  (10) 
shows  that  piezoelectric  crystals  with  high  velocity  of  sound  are  preferable  in  wide-band  acoustic  transducers,  which 
represents  not  an  obvious  conclusion.  The  speed  of  sound  in  piezomateria!  is  not  included  in  consideration  of  standard 
resonantly  detecting  ultrasonic  transducers. 

Table  1  lists  the  most  important  parameters  of  main  piezoelectric  materials  that  could  be  employed  for  wide-band 
ultrasonic  detection.  This  table  shows  that  Lithium  niobate  has  the  highest  value  of  piezoelectric  modulus.  However,  the 
highest  possible  piezoelectric  efficiency  for  wide-band  piezoelectric  detection  can  be  achieved  with  piezoceramics,  such  as 
PZT-5H.  Somewhat  lower  sensitivity  could  be  obtained  with  polyvinyldenefluoride  (PVDF)  or  PVDF  copolymers.  Quartz 
and  lithium  niobate,  materials  widely  used  in  piezoelectric  inductry  have  significantly  lower  sensitivity  for  wide-band 
ultrasonic  detection. 

It  is  important  to  note  that  there  could  be  other  considerations,  which  may  strongly  influence  advisability  of  various 
piezoelectric  materials  for  the  wide-band  detection.  The  specific  electrical  capacity  and  the  acoustic  impedance  of 
piezoelectric  element  are  the  two  most  significant  parameters  that  have  to  be  considered  in  addition  to  the  bandwidth  specific 
piezoelectric  efficiency,  TJac. 

Small  specific  capacity  of  transducer  lower  than  several  pF  will  decrease  sensitivity  of  acoustic  detection,  because 
greater  electrical  capcity  of  electronic  input  circuits  in  preamplifier  will  shunt  the  transducer.  This  factor  is  the  mose 
pronounced  in  detection  of  relatively  low  ultrasonic  frequencies  below  2  MHz.  Therefore,  the  specific  capacity  of 
piezoelement  defines  its  minimal  thickness  and  area.  PZT-5  and  piezoceramic  materails  in  general  possess  the  highest  linear 
capacity,  and  may  be  used  in  microminiature  transducer  designs. 

The  main  disdvantage  of  piezoceramic  materials  for  detection  of  wide-band  ultrasonic  transients  (as  well  as 
piezoelectrics  with  crystal  structure)  is  the  occurrence  of  shear  stress  upon  acoustic  scattering  at  the  boundaries  of 
piezoelectric  element.  Transient  shear  waves  result  in  strong  effect  on  the  transducer  sensitivity  within  the  ultrasonic 
detection  band,  and  thus,  in  deformation  of  the  detected  opto-acoustic  profile.  This  problem  does  not  exist  for  piezoelectric 
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polymers,  such  as  PVDF.  Piezoelectric  composite  materials  also  allow  effective  damping  of  shear  waves  that  might  be 
generated  in  piezoceramics.  Therefore,  piezopolymers  and  PZT-epoxy  composite  materials  certainly  have  noticeable 
advantage  for  use  in  transducers  with  ultrasonic  bandwidth  below  20  MHz.  However,  a  great  deal  of  difficulty  in  fabricating 
piezopolymers  and  composites  with  high  quality  flat  surfaces  makes  them  undesirable  for  transducers  operating  in  ultrasonic 
frequency  range  above  20-30  MHz. 

Table  1.  Properties  of  various  piezoelectric  materials  employed  for  wide-band  ultrasound  detection 


Property  /  Material 

Lithium  Niobate 

Quartz 

PZT-5 

PVDF 

Po>  103  kg/ cm3 

4.63 

2.21 

7.5 

1.8 

C[,  10  3  ml  s 

7.26 

5.57 

4.0 

1.4 

p0q,  10 6  kg /m2s 

33.6 

12.3 

30 

2.5 

g33,  Vm/N 

0.23 

0.058 

0.025 

0.18 

£,  dielectric  permeability 

29 

4.55 

1200 

9.4 

EEq  /  C\ ,  pF  /  mm"' MHz 

0.081 

0.017 

6.1 

0.17 

Q££0g33, 10-9  m3  /  sN 

0.99 

0.75 

27 

3.8 

h^,\09N/C 

5.1 

4.4 

3.0 

0.64 

££q/z33,  109  J  ! nv* 

6.7 

0.78 

95 

0.034 

Therefore,  application  of  easy-to-polish  crystals  of  quartz  and  lithium  niobate  has  significant  technological  merit  in 
the  ultrasonic  frequency  range  above  30  MHz.  Piezoceramic  materials  have  very  high  absorption  in  the  ultrasonic  range 
above  20-30  MHz,  and  thus,  not  used  in  this  ultrasonic  frequency  range.  In  contrast,  piezoceramics  is  a  dominating  material 
for  manufacturing  transducers  operating  in  MHz  and  sub-MHz  ranges,  especially  microminiature  transducers. 


Figure  4.  Acoustic  pressure  generated  by  thermal  noise  of  acoustic  transducer  (operating  in  idler  mode)  as  a  function  of  its 
electric  capacity  presented  for  five  different  ultrasonics  bands  (upper  curve  for  the  bandwidth  of  100  MHz). 
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Similarity  of  acoustic  impedances  of  PVDF  films  and  biological  tissues  makes  PVDF  preferable  for  design  of 
continuously  matching  through  acoustic  tract.  However,  through  matching  is  not  as  important  for  wide-band  ultrasonic 
detection  as  concordance  of  acoustic  impedances  between  piezoelement  and  backing  of  the  transducer.  Furthermore,  acoustic 
pressure  increases  upon  transmission  of  ultrasound  wave  from  medium  with  lower  acoustic  impedance  into  the  medium  with 
higher  acoustic  impedance,  thereby  increasing  sensitivity  of  transducer  in  measuring  opto-acoustic  pressure  generated  in  the 
medium  with  lower  acoustic  impedance  (such  as  biological  tissue).  Low  acoustic  impedance  of  polyvinylidenefluoride 
makes  PVDF  the  material  of  choice  for  a  number  of  applications  in  medical  and  biological  sensing.  Design  of  effective 
backing  for  PVDF  transducers  does  not  represent  a  great  problem  and  makes  detected  opto-acoustic  signals  clear  from 
reverberations  within  the  piezoelement.  For  this  reason,  here  we  made  quantitative  estimates  for  wide-band  ultrasonic 
transducers  made  of  PVDF. 

Minimally  detectable  pressure,  pmi„,  is  presented  in  Fig.  4  as  a  function  of  electric  capacity,  CT,  of  piezoelement 
operating  in  open  circuit  mode.  Taking  into  account  that  the  the  thermal  noise  voltage  is  proportional  to  the  square  root  of 
the  piezoelement  capacity  (see  formula  (7)).  The  electric  capacity  is  inversely  proportional  to  the  thickness  of  piezoelement 
and  the  thickness,  in  turn,  is  defined  from  the  upper  limit  ultrasound  frequency,  then  the  minimally  detectable  pressure  can  be 
found  from  the  following  expression: 
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Fig.  4  demonstrates  that  extremely  low  acoustic  pressure  measured  in  units  of  Pascale  could  be  detected  with 
transducers  operating  in  a  very  wide  frequency  range  of  several  tens  MHz. 


Correct  (undisturbed)  detection  of  opto-acoustic  signals  with  transducers  operating  in  short  circuit  mode  is  possible  only 
within  the  temporal  window  equal  the  time  of  sound  propagation  along  the  thisckness  of  piezoelectric  element.  Electrical 
current  through  the  piezoelement  is  proportional  in  this  window  to  the  vibrational  velicity  of  the  transducer  front  surface, 
while  the  rear  surface  stay  still.  Thus,  detection  of  long  opto-acoustic  signals  requires  substantial  thickness  of  the  transducer 
piezoelemnt.  The  drawback  of  thick  piezoelemnts  is  low  electrical  capacity  and  reduced  sensitivity  of  transducer  resulting 
from  shunting  effect  of  comparable  electrical  capacity  of  preamplifier  and  other  electronic  circutry  of  the  transducer.  On  the 
other  hand,  due  to  the  fact  that  resonance  properties  of  the  piezoelement  (and  reverberations)  do  not  manifest  themselves 
within  this  temporal  window,  the  short  circuit  mode  provides  an  even  sensitivity  characterizticS  in  a  wide  ultrasonic 
frequency  range. 

The  range  of  detectable  ultrasonic  frequencies  is  defined  in  the  short  circuit  mode  by  the  time  of  electric  discharge 
of  piezoelement  capacity  through  the  load  resistor.  The  spectral  reasponse  curve  of  transducer  sensitivity  can  be  presented  in 
this  case  as: 

\$\  =  [\  +  {2nfRCT)2)~1'2  ,  (i2 

The  corresponding  bandwidth  at  the  level  of  0.7  Smax  can  be  expressed  as; 
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As  the  electrical  current  through  the  piezoelement,  Ipe  in  the  short  circuit  mode  is  proportional  to  vibrational 
velocity,  the  expressions  of  Ohm  and  piezoelectric  effect  yield  the  following  voltage  at  the  load  resistor: 

UT  =  Rlpe-RCT833clP'  •  (14) 

The  thermal  noise  is  defined  by  the  load  resistor  and  could  be  found  using  consideration  of  Neiquist: 

CjUj  =  AkToRkf  =  \  AkT0  (15) 

where  A f  =  V3  fupper  is  the  ultrasonic  detection  bandwidth  at  the  level  of  half-maximum. 
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Therefore,  the  voltage  of  thermal  noise,  which  defines  minimally  detectable  pressure, 
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(16) 


is  inversely  proportional  to  the  piezoelement  capacity.  The  expression  (16)  found  from  combining  formulas  (13),  (14),  and 
(15)  closely  resembles  that  for  the  open  circuit  mode  (1 1). 

The  short  circuit  mode  is  justified  for  detection  of  short  opto-acoustic  sisgnals  with  nanosecond  temporal  resolution. 
Crystal  piezomaterials,  such  as  lithium  niobate  and  quartz  are  the  materials  of  choise  for  the  wide-band  transducers  operating 
in  the  short  circuit  mode.  An  important  characteristics  of  these  piezoelectric  crystals  for  transducers  operating  in  high 
ultrasonic  frequency  range  (up  to  several  hundred  MHz)  are  associated  with  low  attenuation  of  acoustic  waves,  negligible 
compared  with  piezoceramic  materials.  Furthermore,  required  flatness  and  parallelizm  of  piezoelement  surfaces  make  crystal 
piezoelectrics  much  preferable  for  high  frequency  acoustic  signals. 

Furthermore,  it  is  possible  to  enhance  effectively  the  sensitivity  of  piezoelectric  transducers  relative  to  the  pressure 
in  media  with  relatively  low  acoustic  impedance,  such  as  biological  tissue.  This  enhancement  can  be  achieved  upon 
transmission  of  ultrasonic  waves  from  a  medium  with  low  acoustic  impedance  into  the  the  medium  with  high  acoustic 
impedance.  Thickness  of  each  layer  is  limited  on  one  hand  by  the  necessity  to  avoid  acoustic  diffraction  (layers  may  not  be 
made  too  thick)  and  on  the  other  hand  by  resonances  of  low  ultrasonic  frequencies  (layers  may  not  be  made  too  thin). 
Nevertheless,  the  pressure  amplitude  could  be  enhanced  3-3.5  times  in  a  system  with  three  layers  with  gradually  increasing 
acoustic  impedance,  so  that  each  transmission  increases  pressure  amplitude  1.5 — 1.8  times  in  the  course  of  acoustic  wave 
propagation  through  each  of  the  two  interfaces. 


4.  ULTIMATE  SENSITIVITY  OF  OPTO-ACOUSTIC  DETECTION 

Let  us  estimate  ultimate  detection  capability  of  wide-band  opto-acoustic  transducers  described  in  Section  3  and 
applied  for  measurements  of  wide-band  pressure  transients  described  in  Section  2.  Here  we  calculate  minimally  detectable 
temperature-rise  and  correspondent  volume  increase  in  a  small  optically  absorbing  sphere  embedded  in  optically  attenuating 
medium  with  acoustic  impedance  matching  that  of  the  sphere.  An  optimal  ultrasonic  detection  band  for  acoustic  waves 
induced  by  laser  pulses  in  this  sphere  is  given  by  expression  (4).  The  acoustic  amplitude  emitted  by  this  sphere  can  be  found 
from  the  expression  (1),  which  could  be  also  presented  as: 

<l7) 
2  r 

where  p0  is  the  densitity  of  the  sphere  and  the  surrounding  medium,  and  AT  is  the  temperature-rise  caused  by  absorption  of 
laser  pulses.  The  pressure  increase  measured  at  the  distance  r  from  the  sphere  is  determined  by  an  instant  pressure  increase 
in  the  heated  sphere  multiplied  by  a  reducing  factor  equal  to  the  ratio  of  the  sphere  radius  and  the  distance  to  the 
measurement  point,  r.  Additional  factor  2  appears  from  the  fact  that  thermoelastic  pressure  is  devided  equaly  between  two 
waves  propagating  outside  from  the  sphere  and  insdie  into  the  sphere.  The  rarefaction  phase  in  the  spherical  acoustic  wave 
could  be  explained  as  a  compression  wave  propagating  into  the  sphere  and  reflecting  from  the  center  with  sign  change. 

For  detection  of  a  spherical  acoustic  wave  with  ultimate  sensitivity,  one  needs  to  employ  a  spherical  acoustic 
transducer  (or  array  transducer  with  equivalent  synthetic  aperture)  shaped  to  match  the  spherical  phase  front  of  the  acoustic 
wave.  The  focus  of  such  a  transducer  needs  to  match  the  location  of  the  absorbing  sphere.  Therefore,  the  area  of  such 
transducer  and  the  radius  of  its  spherical  curvature  could  be  expressed  as: 

A  =  Ur2  (18) 

Bringing  together  the  expressions  (17),  (18),  (4)  and  (9)  one  can  get  the  following  equation  for  minimally  detectable 
temperature  increase  in  laser-irradiated  sphere: 
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where  /l  is  the  piezoelectric  deformation  modulus,  £2  is  the  solid  angle  covered  by  the  acoustic  transducer  and  measured 
from  the  focal  point  located  in  the  center  of  the  absorbing  sphere. 

Formula  (19)  shows  that  minimally  detectable  temperature  is  inversely  proportional  to  the  square  root  of  the  volume 
absorbing  laser  energy.  Thus,  an  ultimate  sensitivity  of  the  opto-acoustic  detection  has  limitation  associated  with  physical 
dimensions  of  detected  object.  Any  increase  in  spatial  resolution  of  the  optoacoustic  imaging  will  unavoidably  result  in  an 
increase  of  minimally  detectable  pressure  or  temperature. 


Figure  5.  Minimally  detectable  temperature  jump  in  tissue  as  a  function  of  characteristic  dimensions  of  absorbing  sphere  for 
three  piezoelectric  materails.  Curve  for  each  piezoelectric  material  is  presented  in  an  optimal  ultrasonic  frequency 
range. 

Figure  5  depicts  minimal  detectable  temperature  generated  in  the  volume  of  a  small  absorbing  sphere  as  a  function 
of  the  sphere  radius  for  wide-band  acoustic  transducers  made  of  various  piezoelectric  materials.  Maximum  sensitivity  could 
be  obtained  with  piezoceramic  transducers.  Lithium  niobate  possesses  somewhat  lower  sensitivity,  and  the  lowest  sensitivity 
can  be  achieved  with  piezopolymers,  such  as  PVDF.  However,  even  with  PVDF  transducers  it  is  possible  to  detect 
impressively  low  temperature  jumps  in  biological  tissues  that  usually  possess  substantial  values  of  thermoelastic  expansion 

coefficient  (2  -10)  *10  4 [1  / AT] ).  Fig.  5  demonstrates  that  with  optimal  design  of  the  spherical  piezoelectric 
transducer  (solid  angle  of  detection  £2  =  1),  it  is  possible  to  detect  a  temperature  increase  of  only  AT- 1-2  pK  in  a  spherical 
tissue  volume  of  about  1-mm3,  which  permits  application  of  low  laser  irradiation  fluences  significantly  below  safe  levels  of 
laser  irradiation  of  human  tissue  [27]. 


5.  DISCUSSION  OF  POSSIBLE  APPLICATIONS 

The  major  trade-off  in  opto-acoustic  imaging  is  between  sensitivity  and  spatial  resolution.  Therefore,  various 
applications  of  the  optoacoustic  tomography  require  different  piezoelectric  transducers  with  ultrasonic  detection  bandwidth 
optimized  for  each  application.  One  of  the  most  prominent  medical  applications  is  the  detection  and  characterization  of  small 
deeply  located  tumors  in  large  human  organs,  such  as  breast  or  brain.  Significant  depth  of  some  tumors  (up  to  6-7  cm  from 
the  surface)  prohibits  propagation  of  high  frequency  ultrasound  to  the  acoustic  transducer.  Therefore,  spatial  resolution  that 
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may  be  achieved  is  about  0.3-mm  to  0.5-mm.  However,  the  emphasis  in  deep  opto-acoustic  imaging  could  be  made  on  high 
sensitivity  in  order  to  detect  small  (2-mm)  tumors  at  the  maximum  possible  depth. 

Assuming  that  a  small  spherical  tumor  is  relatively  transparent  for  optical  radiation,  i.e.  /AaRp  <  1,  and  laser  pulses 

propagate  to  the  tumor  as  spherical  photon  waves  in  tissue  with  effective  optical  attenuation,  JLeff'  one  can  °bta*n  an 
expression  for  the  temperature  increase  inside  spherical  tumor  located  at  a  depth,  Z ,  within  the  irradiated  tissue: 

3/yE(r=  Rp)  exp (-^z) 

P0cp  P0cp  PeffZ 

Substitution  of  (20)  in  (19)  yields  an  equation  that  could  be  employed  in  practical  estimates: 

(Heff 

Let  us  make  our  estimates  using  the  following  parameters,  typical  for  a  small  breast  carcinoma: 

pa  =  0.6  cm'1  =  60  m_1;  £0  =  100  mJ/cm2  =  103  J/m2\  Q  =  l;  (3  =  10~3K~l\  pQcp  =  4*  106  J! m3K 

The  piezoelectric  transducer  made  of  PZT-5  ceramics  has  not  only  high  piezoelectric  pressure  modulus,  but  also 
substantial  acoustic  impedance,  pC/,  which  results  in  two  times  greater  pressure  amplitude  compared  with  the  pressure  in 

2  113 

tissue.  The  combination  of  all  parameters  used  yield  the  following  piezoelectric  efficiency:  ESqH  =  3.8  *  10  Jim  , 
and  therefore 

(PejfZ-max ) exP(^^max )  =  3-2  *  10 6*(R[mm])312  (22) 


^max)exP(^j5f^max)  ~ 


OJp/JLaEo  I^Qh2^R3p 


Tumor  Radius,  mm 


Figure  6.  Maximum  depth  of  a  small  tumor  detection  with  opto-acoustic  imaging  performed  in  forward  mode. 

Maximum  depth  of  opto-acoustic  detection  for  a  small  tumors  in  the  breast  is  shown  in  Fig.  6  as  a  function  of  tumor 
radius.  This  figure  demonstrates  that  a  tumor  with  characteristic  dimension  of  about  1-mm  may  potentially  be  detected  at  the 
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6.  CONCLUSION 
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linearly  with  decreasing  ultrasonic  detection  bandwidth.  The  parameter  which  d'eTnes'omtatl  d"  “pK“n“  an,i 'decreases 
determined  in  this  paper  and  called  rhe  bandwidth  speeilic  pieaoeto  ceffi  tcy  T  <  eT  „  “I"  It  'rh  '  *“ 
includes  the  detection  bandwidth,  the  piezoelectric  modulus  and  the  sn.Td  of  ?  <S"™\(I0»-  Th“  Parameter 
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wide-band  piezoelectric  detection,  we  delermined  three  piezoelectric  material  th  t  °f  P“'“m'ters  'hat  de,in'  sensitivity  of 
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Piezoelectric  Materia  1 


Lithium  Niobate.  LiNbQ3 


Piezoelectric  Ceramics.  PZT-S 
Polyvinyldenefluoride,  PVDF 
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20  kHz  to  20  MHz 


1  MHz  to  80  MHz 
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low  frequency  elements  too  thick 


high  frequency  absorption  too  high 
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